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ABSTRACT 


Illicit  drug-trafficking  is  a  major  concern  of  the  United  States  and  is  a  primary  pillar  of  President 
Barack  Obama’s  Strategy  to  Combat  Transnational  Organized  Crime.  In  the  eastern  Pacific  and 
Caribbean  Sea,  drug-trafficking  organizations  operate  a  variety  of  vessels  to  transit  drugs  from 
South  America  to  the  United  States.  Joint  Interagency  Task  Force  (JIATF)  South,  in  cooperation 
with  partner  agencies  and  nations,  detects,  tracks,  and  interdicts  illegal  drug-trafficking  in  this 
region.  In  this  thesis,  we  develop  a  probability  model  based  on  intelligence  inputs  to  generate  a 
spatial  temporal  heat  map  specifying  the  likely  location  of  targets  over  time.  We  also  formulate  a 
path-finding  model  that  takes  the  heat  map  as  input  and  determines  route  characteristics  through 
transit,  such  as  departure  times,  waypoints,  and  speed.  We  link  the  results  of  our  models  to  a 
separate  effort  that  seeks  to  provide  JIATF  South  with  an  optimal  search  plan  to  maximize  the 
expected  amount  of  drugs  seized.  We  show  that  our  path-finding  model  accurately  reproduces 
the  target’s  track  even  when  the  target  transits  along  a  complex  route.  Furthermore,  we  show 
that  the  optimal  search  plan  based  on  our  path-finding  model  is  nearly  identical  to  the  search 
plan  using  the  known  parameters. 
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Executive  Summary 


Illicit  drug-trafficking  originating  from  South  America  is  a  major  concern  of  the  United  States, 
and  continued  expansion  of  drug-trafficking  significantly  affects  the  national  security  interests 
of  the  United  States.  Organizations  tasked  with  combatting  drug  traffickers  are  faced  with 
several  problems.  They  must  detect,  track,  and  interdict  illegal  drug  traffickers  en  route  to 
various  follow-on  transit  zones  in  Central  America,  Mexico,  and  the  Caribbean.  They  utilize 
incoming  intelligence  from  various  sources  to  determine  where  and  what  assets  to  employ. 
Their  efforts,  while  productive,  are  not  able  to  stop  a  large  majority  of  the  drugs  from  reaching 
the  United  States.  They  battle  Mexican  and  Colombian  Drug  Trafficking  Organizations  (DTOs) 
spread  throughout  Latin  America.  The  DTOs  have  extensive  capabilities  and  access  to  large 
amounts  of  money.  Combined  with  the  close  proximity  to  the  United  States,  the  DTOs  pose  an 
extreme  threat  to  the  nation’s  security. 

Operators  at  Joint  Interagency  Task  Force  (JIATF)  South  obtain  various  types  of  intelligence 
information  and  utilize  that  information  when  creating  search  and  interdiction  plans.  They  may 
have  at  their  disposal  information  concerning  departure  times,  departure  and  arrival  locations, 
type  of  vessel,  speed,  drug  load,  and  waypoints.  At  all  times,  there  is  a  level  of  uncertainty  as¬ 
sociated  with  the  intelligence,  which  can  be  represented  by  a  probability  distribution  for  each  of 
the  aforementioned  pieces  of  information.  Given  this  information  and  surrounding  uncertainty, 
JIATF  South  operators  can  benefit  from  a  visual  heat  map  displaying  the  likely  route  of  a  target 
or  targets  over  time.  In  this  thesis,  we  address  the  problem  of  creating  a  heat  map  of  likely  drug 
trafficker  locations  based  on  intelligence  inputs. 

One  approach  to  generating  a  heat  map  is  through  simulation.  Simulations  have  the  benefit  of 
being  able  to  account  for  many  real  world  complexities  like  weather.  However,  they  also  may 
be  computationally  expensive  and  difficult  to  validate  and  verify.  We  desire  to  complement 
and  vet  such  complicated  simulation  by  developing  more  analytically  tractable  models.  We 
develop  probability  models  to  generate  a  heat  map  to  provide  operators  with  better  situational 
awareness  about  the  likely  locations  of  targets,  both  in  the  current  period  and  future  periods. 
We  first  build  a  simple  probability  model  that  accounts  for  uncertainty  only  in  the  departure 
and  arrival  locations.  Using  this  model  as  a  baseline,  we  further  enhance  its  functionality  to 
account  for  uncertain  speed  and  time  of  departure,  as  well  as  consideration  for  multiple  legs, 
multiple  targets  in  the  same  region,  and  the  ability  to  output  expected  amount  of  drugs.  Given 
the  probability  distributions  for  the  intelligence  inputs  (e.g.,  speeds,  locations,  and  times),  the 


model  produces  the  spatial  temporal  probability  map  of  the  trafficker’s  path.  Using  the  heat 
map  as  input,  we  then  address  the  problem  of  determining  a  drug  trafficker’s  transit  route  based 
on  the  collected  intelligence  information. 

After  developing  a  model  to  transform  intelligence  inputs  into  a  heat  map,  we  next  turn  to  the 
inverse  problem:  Given  a  heat  map,  what  are  the  most  likely  intelligence  inputs  used  to  create  it. 
We  adopt  this  approach  to  compare  the  path-finding  model’s  results  with  some  known  inputs  via 
the  search  vehicle  optimization  model.  Additionally,  it  is  possible  that  the  only  data  available 
for  JIATF  South  planners  will  be  a  heat  map,  thus  providing  a  means  to  estimate  intelligence 
parameters  is  beneficial  for  the  planners. 

We  produce  the  necessary  intelligence  parameters  through  the  use  of  various  regression  tech¬ 
niques  and  consideration  for  the  input  requirements  of  the  optimization  model.  The  path-finding 
model  approximates  the  path  taken  by  the  target  as  a  sequence  of  piecewise  linear  legs,  where 
the  collection  of  all  legs  produces  the  target’s  entire  track  from  departure  location  to  arrival 
location.  Standard  regression  techniques  are  incapable  of  automatically  identifying  where  a 
new  leg  begins,  which  occurs  when  speed  or  heading  changes.  As  a  result,  we  employ  more 
advanced  regression  techniques,  namely  Multivariate  Adaptive  Regression  Splines  (MARS).  In 
accordance  with  the  search  vehicle  optimization  model,  we  define  a  leg  as  a  segment  of  the 
target’s  track  that  has  a  constant  heading  and  speed.  This  creates  a  piecewise  linear  path  for  the 
target.  Specifically,  we  determine  the  target’s  departure  time  range,  expected  departure  time, 
arrival  time  range,  expected  arrival  time,  speed,  departure  and  arrival  locations,  and  waypoint 
locations. 

The  final  aspect  of  the  path-finding  model  is  specifying  the  Area  of  Uncertainty  (AOU)  of  a  tar¬ 
get  as  it  traverses  along  a  path  leg.  Roughly  speaking,  the  AOU  is  the  size  of  the  target’s  "blob" 
on  the  heat  map.  The  AOU  along  a  target’s  track  represents  the  uncertainty  in  both  departure 
speed  and  departure  time.  With  these  uncertainties,  an  area  exists  in  which  it  is  possible  for  the 
target  to  be  present  at  any  point  in  the  associated  area.  The  search  plan  algorithm  requires  that 
the  AOU  be  constant  on  each  leg  of  the  path.  Thus,  if  the  AOU  changes  significantly,  we  break 
one  leg  of  the  path  into  multiple  legs  even  if  the  speed  and  heading  do  not  change.  We  develop 
some  heuristics  to  accomplish  this. 

We  develop  two  case  studies  and  the  corresponding  intelligence  information  from  which  we 
generate  a  heat  map  according  to  the  probability  model.  Next,  we  use  our  path-finding  algo¬ 
rithm  to  estimate  the  input  parameters  and  compare  how  well  they  match  to  the  true  values. 


xviii 


We  show  the  capability  to  identify  speed  changes,  changes  in  heading  or  course,  and  changes 
in  the  AOU  surrounding  the  target  throughout  the  path.  We  illustrate  the  models’  functionality 
with  two  case  studies  of  varying  complexity.  First,  a  simple  straight  path  between  the  origin 
and  destination  shows  the  path-finding  model’s  ability  to  identify  the  target’s  starting  location, 
destination,  speed,  and  associated  travel  times.  Combining  these  route  characteristics,  we  are 
able  to  reproduce  the  target’s  likely  track  from  start  to  finish.  Supplying  our  path-finding  results 
into  the  search  vehicle  allocation  model  shows  similar  search  plan  results  when  compared  with 
the  original  case  parameters  used  in  the  probability  model.  Next,  we  demonstrate  a  more  com¬ 
plex  scenario  with  a  target  traveling  to  the  destination  with  a  course  change  during  the  journey. 
Our  model  accurately  identifies  this  course  change  and  reproduces  the  target’s  track.  These  two 
case  studies  are  representative  of  actual  drug  trafficker  routes  but  additional  analysis  is  needed 
to  validate  the  model’s  ability.  Finally,  we  provide  the  search  vehicle  allocation  model  with  the 
path-finding  model  results  for  a  scenario  with  two  targets.  The  resulting  search  plan  is  nearly 
identical  to  the  one  generated  when  the  true  intelligence  inputs  are  known.  This  provides  strong 
evidence  that  we  have  provided  a  key  link  between  the  probability,  path- finding,  and  optimiza¬ 
tion  models. 
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CHAPTER  1: 
Introduction 


Illicit  drug-trafficking  is  a  major  concern  of  the  United  States  and  is  a  primary  pillar  of  Pres¬ 
ident  Barack  Obama’s  Strategy  to  Combat  Transnational  Organized  Crime.  Transnational  Or¬ 
ganized  Crime  (TOC)  encompasses  many  areas  of  global  crime  to  include  human  smuggling, 
trafficking  in  persons,  weapons  trafficking,  intellectual  property  theft,  cybercrime,  and  drug¬ 
trafficking  (National  Security  Staff,  2011).  Drug-trafficking  has  expanded  in  recent  years,  and 
organizations  have  become  increasingly  more  violent  and  less  centralized.  The  U.S.  govern¬ 
ment  lists  Drug  Trafficking  Organizations  (DTOs)  as  the  "greatest  organized  crime  threat"  to 
the  United  States  (National  Drug  Intelligence  Council  [NDIC],  2008).  Murder  rates  in  known 
drug-trafficking  countries  have  increased,  and  continued  expansion  of  drug-trafficking  signif¬ 
icantly  affects  the  national  security  interests  of  the  United  States.  This  thesis  contributes  to 
the  ongoing  research  to  improve  the  effectiveness  of  the  agencies  tasked  with  combatting  il¬ 
licit  drug-trafficking.  Based  on  incoming  intelligence,  the  analysis  models  movements  of  drug 
traffickers. 


1.1  Background 

Latin  America  is  a  prominent  player  in  the  global  drug  market,  with  involvement  in  nearly 
all  areas  of  the  illicit  drug  supply  chain.  Activities  in  the  region  include  drug  crop  farming, 
drug  production,  drug-trafficking,  and  drug  consumption.  The  South  American  countries  in  the 
Andean  region,  particularly  Colombia,  Peru,  and  Bolivia,  are  the  sole  producers  of  cocaine  for 
the  global  market  (United  States  Department  of  State,  2012).  United  States  Department  of  State 
(2012)  estimates  total  cocoa  production  in  these  three  countries  at  187,500  hectares  in  2010. 
This  is  down  from  a  peak  in  2007  of  232,500  hectares.  Combining  the  Andean  region’s  ability 
to  produce  cocaine  and  their  close  proximity  to  the  United  States,  makes  Latin  America  the 
central  hub  of  the  illicit  drug  market  (Seelke,  Wyler,  Beittel,  &  Sullivan,  2011).  Additional 
factors  that  contribute  to  the  region’s  role  in  the  drug  market  include  poverty,  inequality,  the 
absence  of  other  fiscal  opportunities,  and  a  weak  governing  force  (Seelke  et  al.,  2011).  These 
factors  have  left  DTOs  nearly  free  to  act  as  they  wish  in  the  supply  of  illicit  drugs. 

A  variety  of  DTOs  operate  in  Latin  America,  and  they  range  from  small  family-based  operations 
to  global  networks  that  are  well  connected  throughout  the  world  to  include  ties  with  terrorist 
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organizations.  The  more  advanced  DTOs  employ  considerable  paramilitary  and  counterintelli¬ 
gence  assets  to  contend  with  legitimate  police  forces  and  conduct  operations  with  minimal  risk 
via  corrupt  governments  (Seelke  et  al.,  201 1).  The  family-based  operations  typically  specialize 
in  a  few  areas  of  the  drug  supply  chain. 

Mexican  and  Colombian  DTOs  are  the  most  prominent  in  the  region.  They  currently  have 
operations  worldwide  to  include  Europe,  West  Africa,  and  the  Asia-Pacific  regions  as  well  as 
connections  to  organized  crime  and  terrorists  (Seelke  et  al.,  2011).  The  NDIC  (2008)  reports 
that  Mexican  and  Colombaian  DTOs  generate,  remove,  and  launder  between  $18  billion  and  $39 
billion  in  wholesale  drug  proceeds  annually.  Currently,  Mexican  DTOs  are  the  dominant  groups 
supplying  most  illicit  drugs  to  the  United  States.  Seven  Mexican  DTOs  are  currently  operating 
in  all  nine  of  the  Organized  Crime  Drug  Enforcement  Task  Forces  (OCDETF)  regions.  Table  1 . 1 
and  Figure  1.1  show  the  Mexican  DTOs  and  their  primary  drugs  and  regions  of  operations. 


Figure  1.1:  Nine  regions  of  the  Organized  Crime  Drug  Enforcement  Task  Force.  All  seven 
Mexican-based  DTOs  operate  in  all  of  the  OCDETF  regions.  (From  NDIC,  2011) 


Over  the  last  ten  years,  the  Mexican  DTOs  have  replaced  many  of  the  Colombian  DTOs  that 
controlled  the  market  in  the  1980s  and  1990s.  The  Mexican  organizations  have  a  geographical 
advantage  over  the  Colombian  drug  traffickers.  They  have  access  and  control  of  the  smuggling 


2 


DTO 

Primary  Drugs 

Primary  Regions 

Sinaloa  Cartel 

Cocaine 

Heroin 

Marijuana 

MDMA 

Methamphetamine 

Florida/Caribbean 
Great  Fakes 

Mid-Atlantic 

New  England 

New  York/New  Jersey 
Pacific 
Southeast 
Southwest 

West  Central 

Los  Zetas 

Cocain 

Marijuana 

Florida/Caribbean 
Great  Fakes 
Southeast 
Southwest 

Gulf  Cartel 

Cocaine 

Marijuana 

Florida/Caribbean 

Mid-Atlantic 

New  England 

New  York/New  Jersey 
Southeast 
Southwest 

Juarez  Cartel 

Cocaine 

Marijuana 

Great  Fakes 

New  York/New  Jersey 
Pacific 
Southeast 
Southwest 

West  Central 

BLO 

Cocaine 

Heroin 

Marijuana 

Southeast 

Southwest 

LFM 

Cocaine 

Heroin 

Marijuana 

Methamphetamine 

Southeast 

Southwest 

Tijuana  Cartel 

Cocaine 

Heroin 

Marijuana 

Methamphetamine 

Great  Fakes 

Pacific 

Southwest 

Table  1.1:  Mexican-based  DTOs  operate  in  all  nine  of  the  OCDETF  regions  and  supply  most 
illicit  drugs  to  the  United  States.  (From  NDIC,  201 1) 
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routes  along  the  U.S.  southwest  border,  and  these  routes  or  plazas  are  fiercely  contested  among 
the  Mexican  DTOs  (NDIC,  2011).  Additionally,  the  ability  to  obtain,  transport,  and  distribute 
the  illicit  drugs  in  demand  in  the  United  States  allows  them  to  dominate  the  market.  The  Sinaloa 
Cartel,  with  operations  in  all  nine  OCDETF  regions,  is  the  principal  Mexican  DTO.  It  boasts 
an  extensive  transportation  and  distribution  network,  a  vast  network  of  personnel  to  aid  in  traf¬ 
ficking  operations,  and  control  of  heroin,  marijuana,  and  methamphetamine  production  (NDIC, 
2011).  According  to  NDIC  (2011),  the  Mexican  DTOs  will  remain  the  dominant  U.S.  illicit 
drug  suppliers  for  the  near  future. 

The  Colombian  DTOs,  the  Cali  cartel  and  Medellin  cartel  that  dominated  the  trade  in  the  1980s 
and  1990s,  have  been  replaced  by  Norte  del  Valle  DTO.  Rather  than  being  one  central  cartel 
such  as  their  predecessors,  the  Norte  del  Valle  consists  of  many  smaller  cells  and  is  less  hi¬ 
erarchical  and  more  decentralized  (Seelke  et  al.,  2011).  State  Department  designated  Foreign 
Terrorist  Organizations  (FTOs),  the  Revolutionary  Armed  Forces  of  Colombia  (FARC),  the 
National  Fiberation  Army  (EFN),  and  the  United  Defense  Forces  of  Colombia  (AUC),  control 
and  finance  Colombian  illicit  drug-trafficking  operations.  The  FTOs  then  use  the  drug  money 
to  attack  Colombian  police  and  military  (United  States  Department  of  State,  2012).  These  or¬ 
ganizations  are  becoming  more  involved  in  drug-trafficking  operations.  They  tax  cocoa  farmers 
and  buyers  and  have  controlling  interests  in  various  operations  of  the  drug  supply  chain  (Seelke 
et  al.,  2011).  When  the  AUC  demobilized,  smaller  paramilitary  groups  formed  and  increased 
their  involvement  in  the  drug  market,  proving  to  be  a  major  challenge  for  law  enforcement. 
Their  increased  presence  has  caused  a  rise  in  violence  in  both  Colombian  cities  and  rural  vil¬ 
lages  (United  States  Department  of  State,  2012). 

Many  of  the  drugs,  specifically  cocaine,  that  reach  U.S.  soil  transit  from  South  America  through 
the  Central  America-Mexico  transit  routes  or  via  the  Caribbean  transit  routes  en  route  to  the 
United  States.  The  Consolidated  Counterdrug  Database  (CCDB)  documented  1,494,  1,228,  and 
804  metric  tons  of  cocaine  movement  through  the  Western  Hemisphere  Transit  Zone  in  fiscal 
years  2008,  2009,  and  2010,  respectively  (Office  of  National  Drug  Control  Policy,  2013).  The 
primary  transportation  route  is  through  the  Central  America-Mexico  corridor  to  include  terri¬ 
torial  waters  in  both  the  eastern  Pacific  and  the  Caribbean  Sea.  This  drug  flow  accounts  for 
95%  of  the  cocaine  going  to  the  United  States.  More  than  half  of  these  drugs  transit  through 
Central  America  first  (Seelke  et  al.,  2011).  This  is  a  shift  from  the  operations  of  the  1980s 
and  early  1990s,  when  drugs  were  mainly  smuggled  through  the  Caribbean  to  Florida  (Seelke 
et  al.,  2011).  In  2010,  5%  of  the  cocaine  transiting  to  the  United  States  was  via  the  Caribbean, 
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whereas  in  2001,  26%  was  transiting  via  the  Caribbean  (United  States  Senate  Caucus  on  Inter¬ 
national  Narcotics  Control,  2012).  In  September  2011,  President  Obama  identified  22  countries 
as  major  drug  transit  or  major  illicit  drug  producing  countries.  Of  the  22  identified,  19  of  the 
countries  are  in  Latin  America  or  the  Caribbean  (United  States  Department  of  State,  2012).  In 
the  Caribbean  region  specifically,  The  Bahamas,  the  Dominican  Republic,  Haiti,  and  Jamaica 
are  major  countries  used  in  the  illegal  trafficking  of  cocaine  (United  States  Senate  Caucus  on  In¬ 
ternational  Narcotics  Control,  2012).  Approximately  one  half  of  the  cocaine  that  flows  through 
the  Caribbean  transits  through  the  Dominican  Republic  (United  States  Senate  Caucus  on  In¬ 
ternational  Narcotics  Control,  2012).  Puerto  Rico  presents  another  major  transit  hub.  Since 
the  country  is  within  the  U.S.  customs  zone,  regulations  to  get  into  the  United  States  are  less 
stringent,  thus  70%  to  80%  of  the  cocaine  moving  through  Puerto  Rico  is  bound  for  the  United 
States  while  the  remaining  cocaine  is  consumed  within  Puerto  Rico  (United  States  Senate  Cau¬ 
cus  on  International  Narcotics  Control,  2012).  The  Caribbean  presents  three  main  corridors  in 
Figure  1.2  for  the  illicit  transport  of  drugs. 


Figure  1.2:  Illicit  Caribbean  drug  routes  are  divided  into  three  main  corridors:  Central,  Eastern, 
and  ABC.  (From  United  States  Senate  Caucus  on  International  Narcotics  Control,  2012) 

The  central  corridor  is  a  primary  avenue  for  moving  cocaine  from  South  America  through  Ja¬ 
maica,  Haiti,  the  Dominican  Republic,  or  The  Bahamas  into  the  United  States.  The  eastern 
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corridor  moves  heroin  and  cocaine  through  eastern  Caribbean  nations,  Trinidad,  and  Tobago 
to  Puerto  Rico  and  ultimately  into  the  United  States.  The  ABC  corridor  of  Aruba,  Bonaire, 
and  Curacao  is  prevalent  among  cocaine  traffickers  as  well  (United  States  Senate  Caucus  on 
International  Narcotics  Control,  2012). 

In  2011,  80%  of  drugs  destined  for  the  United  States  moved  via  the  waterways  and  the  remain¬ 
ing  20%  via  aircraft  (United  States  Senate  Caucus  on  International  Narcotics  Control,  2012). 
Maritime  and  air  transit  routes  from  Joint  Interagency  Task  Force  (JIATF)  South  are  shown 
in  Figures  1.3  and  1.4.  DTOs  using  maritime  Central  America-Mexico  transit  routes  utilize 
both  Caribbean  and  eastern  Pacific  routes  as  evident  in  Figure  1.3.  The  air  routes  leaving  the 
Colombia- Venezuela  border  avoid  Colombian  air  space  in  an  effort  to  evade  surveillance  by 
the  U.S.  Government  and  Colombian  Government  as  part  of  the  Air  Bridge  Denial  Agreement. 
Colombian  officials  cite  that  the  surveillance  and  interdiction  efforts  force  smugglers  to  remain 
clear  of  Colombian  air  space  (GAO,  2009). 


Figure  1.3:  Maritime  drug  transit  routes  in  2011.  80%  of  the  drugs  destined  for  the  United 
States  transit  via  waterways.  (From  United  States  Senate  Caucus  on  International  Narcotics 
Control,  2012) 


The  primary  means  of  transport  is  via  maritime  vessels,  which  include  go-fast  boats,  pleasure 
craft,  fishing  vessels,  and  Self-propelled  Semi-submersible  (SPSS).  The  Caribbean  has  seen  a 
significant  increase  in  the  number  of  go-fast  boats  used  to  transport  illicit  drugs.  These  vessels 
travel  at  high  speeds,  are  small,  display  nearly  no  radar  cross-section,  and  are  hard  to  see 
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Figure  1 .4:  Air  drug  transit  routes  in  201 1 .  (From  United  States  Senate  Caucus  on  International 
Narcotics  Control,  2012) 


during  the  day,  all  of  which  makes  them  difficult  to  find,  track,  and  interdict.  They  can  carry 
as  much  as  two  metric  tons  of  cocaine  and  their  transits  usually  last  one  or  two  days  (Office  of 
National  Drug  Control  Policy,  2013).  In  the  eastern  Pacific,  fishing  vessels  typically  depart  from 
Colombian  and  Ecuadorian  ports  carrying  several  tons  of  cocaine  to  ports  in  Central  America 
or  Mexico.  These  operations  can  last  up  to  six  weeks  and  inject  substantial  amounts  of  cocaine 
into  the  Central  America-Mexico  transit  route.  Drug  traffickers  have  recently  been  employing, 
to  a  greater  extent,  more  subtle  maritime  transportation  in  the  form  of  SPSSs. 

The  use  of  SPSSs,  such  as  the  one  in  Figure  1.5,  to  transport  cocaine  is  on  the  rise  in  the 
Caribbean  and  eastern  Pacific.  Between  2001  and  2007,  23  events  involving  SPSSs  were  doc¬ 
umented,  compared  to  62  SPSS  events  in  the  first  three  quarters  of  fiscal  year  2008  (USCG, 
2009).  SPSSs  are  built  mainly  from  wood  and  fiberglass  in  the  jungles  of  Colombia.  They 
are  capable  of  carrying  up  to  10  metric  tons  of  cocaine  with  a  range  of  2,500  nautical  miles. 
They  range  in  length  from  25  to  80  feet  and  travel  up  to  13  knots.  A  crew  of  4  to  5  operates 
the  vessel,  which  is  designed  to  scuttle  upon  detection  by  law  enforcement  (USCG,  2009).  Re¬ 
cently,  DTOs  have  been  improving  their  capabilities  with  the  construction  of  fully-submersible 
vessels.  In  July  of  2010,  the  Ecuador  police  seized  a  drug  submarine  while  under  construction 
in  the  jungles  of  Ecuador  that  is  100  feet  long  and  capable  of  diving  to  65  feet  (Kraul,  2010). 
The  switch  to  fully  submersible  submarines  poses  a  significantly  harder  threat  to  combat. 
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Figure  1.5:  SPSS  as  seen  by  U.S.  Coast  Guard  prior  to  being  boarded.  (From  USCG,  2009) 


Antidrug  efforts  have  been  a  mainstay  of  U.S.  drug  policy  in  Latin  America  and  the  Caribbean 
for  the  last  30  years,  with  the  United  States  spending  $14.6  billion  to  fight  illegal  drugs.  The 
strategy  is  to  stop  it  at  the  source  by  destroying  drug  producing  crops,  interrupting  shipments, 
and  improving  socio-economic  factors  in  countries  prone  to  DTOs.  The  programs  are  funded  by 
the  U.S.  Department  of  State,  U.S.  Agency  for  International  Development,  and  the  Department 
of  Defense  (DoD).  The  DoD  works  with  foreign  countries  to  train,  equip,  and  improve  counter 
drug  efforts.  National  Security  Directive  18,  signed  by  President  Bush  on  August  21,  1989,  and 
the  Andean  Initiative  put  counter  narcotics  on  the  DoD’s  list  of  core  priorities.  It  is  now  the  lead 
federal  agency  responsible  for  detection  and  monitoring  of  illegal  drugs  from  foreign  countries 
into  the  United  States.  The  National  Defense  Authorization  Acts  of  1991  and  1998  as  amended 
through  FY2011  allows  the  DoD  to  provide  support  to  certain  countries,  mainly  Colombia. 
More  recently,  the  Andean  Counterdrug  Initiative  (ACI)  proposed  by  President  Clinton  in  2000 
and  passed  through  Congress  increased  funding  for  counter  drug  operations  in  South  America. 
The  primary  countries  of  the  ACI  include  Colombia,  Peru,  and  Bolivia  (Seelke  et  al.,  2011).  In 
fiscal  year  2008,  the  Merida  Initiative  increased  funding  for  counter  drug  operations  in  Mex¬ 
ico  and  since  then  efforts  have  shifted  from  Colombia  to  Mexico,  Central  America,  and  the 
Caribbean.  Two  recent  initiatives,  the  Central  America  Regional  Security  Initiative  (CARSI) 
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and  the  Caribbean  Basin  Security  Initiative  (CBSI)  have  brought  more  focus  to  the  region  and 
DoD  operations  specifically  support  these  initiatives  (Seelke  et  al.,  2011). 

1.2  U.S.  Southern  Command 

DoD  counternarcotics  policy  is  established  by  the  Office  of  the  Deputy  Assistant  Secretary 
of  Defense  for  Countemarcotics  and  Global  Threats.  Two  combatant  commands  in  the  west¬ 
ern  hemisphere  execute  the  mission  of  providing  counter  drug  support,  U.S.  Southern  Com¬ 
mand  (SOUTHCOM)  and  U.S.  Northern  Command  (NORTHCOM).  NORTHCOM’s  Area  of 
Responsibility  (AOR)  covers  Mexico  and  SOUTHCOM’s  AOR  covers  all  of  Latin  America 
south  of  Mexico.  SOUTHCOM’s  efforts  to  counter  illegal  drug-trafficking  directly  supports 
the  first  objective  of  their  Command  Strategy  2020,  defend  the  United  States  (Fraser,  2010). 
SOUTHCOM  is  composed  of  five  component  commands,  three  joint  task  forces,  and  one  direct 
reporting  unit.  SOUTHCOM  collaborates  with  partner  nations  to  combat  TOC  in  their  AOR. 
Illicit  drug-trafficking  is  their  primary  focus  as  this  is  a  means  for  TOCs  to  strengthen  their 
influence  in  the  region.  SOUTHCOM  provides  support  to  partner  nations  enabling  them  to  pro¬ 
vide  better  security  in  the  region.  This  support  is  directed  to  interrupt  the  flow  of  illegal  drugs, 
remove  TOC  networks,  and  provide  assistance  to  remove  criminal  influences  in  lawless  areas. 
As  counter  narcotics  are  a  core  priority  of  the  DoD,  SOUTHCOM  is  the  lead  agency  in  the 
AOR  in  carrying  out  this  mission.  Their  efforts  are  in  line  with  the  objectives  of  the  Strategy 
to  Combat  Transnational  Organized  Crime,  the  CARSI,  and  the  CBSI  (United  States  Southern 
Command,  2013). 

1.3  Joint  Interagency  Task  Force  South 

JIATF  South  is  a  task  force  of  SOUTHCOM  that  directly  performs  SOUTHCOM’s  mission  to 
combat  illicit  drug-trafficking.  Located  in  Key  West,  FL,  JIATF  South  detects,  tracks,  and  inter¬ 
dicts  illegal  drug-trafficking  activity  in  the  eastern  Pacific,  Gulf  of  Mexico,  and  the  Caribbean. 
They  work  with  partner  nations,  military  organizations,  and  U.S.  law  enforcement  agencies 
in  carrying  out  their  mission.  U.S.  military  and  partner  nations’  maritime  and  air  assets  are 
employed  in  the  detection,  monitoring,  and  interdiction  of  illicit  drug-trafficking  in  the  region 
(United  States  Southern  Command,  2013).  The  assets  available  to  JIATF  South  can  vary  over 
time  and  have  recently  been  in  decline.  Typically,  U.S.  Navy,  U.S.  Coast  Guard,  and  partner 
nation  maritime  vessels  patrol  the  Caribbean  Sea,  Gulf  of  Mexico,  and  the  eastern  Pacific.  A 
U.S.  Coast  Guard  Law  Enforcement  Detachment  embarks  on  a  U.S.  vessel  to  take  control  dur¬ 
ing  boarding  operations  to  seize  illicit  drugs  and  apprehend  suspects  (United  States  Southern 
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Command,  2013).  In  the  air,  U.S.  military,  interagency,  and  partner  nation  aircraft  execute  de¬ 
tection  and  monitoring  missions  and  direct  interdiction  assets  to  suspicious  surface  vessels  for 
further  investigation. 

In  2011,  JIATF  South  interdicted  119  metric  tons  of  cocaine,  arrested  355  drug  traffickers, 
and  seized  70  aircraft  and  vessels,  which  resulted  in  a  $7.1  billion  loss  of  revenue  for  DTOs 
(Posture  Statement  of  General  Douglas  M.  Fraser,  2012).  They  projected  775  -  930  metric 
tons  of  cocaine  to  move  towards  the  United  States  in  2012  and  in  January,  JIATF  South  be¬ 
gan  Operation  Martillo  with  the  goal  of  disrupting  drug-trafficking  along  Central  American 
coastal  waters  (Posture  Statement  of  General  Douglas  M.  Fraser,  2012).  U.S.  Naval  Forces 
Southern  Command,  in  Mayport,  FL,  supported  Operation  Martillo  with  seven  frigates,  one  re¬ 
plenishment  ship,  and  four  fixed-wing  aviation  squadrons  (Posture  Statement  of  General  John 
F.  Kelly,  2013).  In  the  first  year  of  Operation  Martillo,  JIATF  South  and  partner  nations  have 
disrupted  152  metric  tons  of  cocaine,  confiscated  $7.2  million  in  bulk  cash,  and  seized  4  aircraft, 
8  SPSSs,  and  89  vessels  (Posture  Statement  of  General  John  F.  Kelly,  2013).  This  success  is 
closely  dependent  on  support  from  partner  nations  who  played  a  role  in  67%  of  the  disruptions 
associated  with  Operation  Martillo  (Posture  Statement  of  General  John  F.  Kelly,  2013).  JIATF 
South  has  played  a  major  role  in  the  planing  and  execution  of  this  operation.  Based  on  the 
amount  interdicted  in  2011,  approximately  85%  of  the  cocaine  is  not  stopped.  They  oversee  a 
large  operational  area  and  effective  and  efficient  employment  of  available  assets  is  paramount 
to  their  success. 

1.4  Mathematical  Modeling  Approach 

This  thesis  addresses  some  of  the  uncertainty  involved  with  detection  and  monitoring  of  illicit 
drug-trafficking  in  the  Caribbean  and  eastern  Pacific.  JIATF  South  is  the  primary  agency  that 
is  most  interested  in  this  work.  Given  the  number  of  assets  they  have  and  the  vast  amount 
of  area  under  their  watch,  their  success  is  heavily  dependent  on  employing  available  assets  as 
effectively  as  possible.  The  main  purpose  of  this  overall  endeavor  is  to  provide  better  situational 
awareness  to  operators,  analysts,  and  targeteers  about  the  likely  location  of  targets  in  the  present 
and  future.  To  meet  this  goal,  we  formulate  two  models. 

JIATF  South  faces  two  primary  problems  in  this  context.  First,  the  locations  of  the  targets 
are  unknown.  A  probability  model  creates  a  heat  map  from  various  intelligence  inputs,  which 
include  departure  location,  arrival  location,  and  departure  time.  This  provides  a  better  picture 
of  the  overall  AOR.  Next,  applying  a  path-finding  algorithm  provides  JIATF  South  with  a 
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probable  track  of  the  target,  which  will  aid  in  the  asset  allocation  process.  Employment  of  their 
assets  represents  the  second  problem.  After  developing  an  estimate  for  the  likely  location  of 
a  target  over  time,  where  to  send  their  assets  is  the  next  complexity  to  overcome.  Pietz  and 
Royset  (2013)  develop  an  optimization  model  that  takes  as  input  the  output  from  the  path¬ 
finding  algorithm  and  creates  a  search  plan  for  JIATF  South. 

Based  on  incoming  intelligence,  the  probability  model  creates  a  time-layered  map  indicating 
the  probability  associated  with  a  drug  trafficker  being  in  a  specific  area  at  a  specific  time.  As 
additional  intelligence  is  received,  the  map  is  updated  with  the  new  information.  Following  the 
creation  of  the  probability  map,  an  algorithm  determines  the  target’s  course,  speed,  waypoints 
(if  applicable),  time  of  departure,  and  time  of  arrival.  The  optimization  model  uses  the  algorithm 
output  to  maximize  the  effectiveness  of  JIATF  South  assets  with  respect  to  the  potential  amount 
of  drugs  seized.  With  the  connection  of  our  probability  model  and  path-finding  algorithm  to  the 
optimization  model,  JIATF  South  receives  an  all-inclusive  tactical  decision  aid  to  aid  in  their 
efforts  to  counter  illicit  drug-trafficking. 

1.5  Thesis  Structure 

This  thesis  develops  a  method  to  determine  a  drug  trafficker’s  course  and  speed  through  time 
given  some  initial  intelligence  reports  and  links  additional  work  being  performed  in  support 
of  JIATF  South’s  mission.  Chapter  2  details  a  literature  review  to  support  the  efforts  of  this 
study.  Chapter  3  describes  the  probability  model,  key  uncertainties,  and  results  of  the  probabil¬ 
ity  model.  Chapter  4  details  the  target  path-finding  model  and  key  inputs  to  the  search  vehicle 
optimization  model.  Chapter  5  presents  the  conclusions  and  recommendations  of  this  thesis. 
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CHAPTER  2: 
Literature  Review 


This  thesis  covers  a  range  of  topics  that  have  been  extensively  researched  in  the  past  and  lever¬ 
ages  techniques  from  probability,  regression,  and  optimization.  Topics  such  as  search  theory 
can  be  traced  back  to  World  War  II  and  the  work  of  Koopman  and  colleagues  on  the  Operations 
Evaluation  Group  (OEG).  More  recently,  researchers  at  the  Naval  Research  Laboratory  (NRL) 
have  developed  a  probability  model  of  likely  pirate  locations  based  on  both  intelligence  and 
historical  weather  data.  A  major  contribution  of  this  thesis  links  our  path-finding  model  to  an 
JIATF  South  asset  optimization  model.  This  research  and  other  applicable  work  is  reviewed  to 
give  context  to  this  thesis. 

2.1  Search  Theory 

The  discipline  of  search  theory  originates  from  the  need  to  locate  enemy  submarines  in  World 
War  II.  Koopman  (1946)  presents  initial  findings  in  this  area  of  research  focusing  on  various 
types  of  optimal  search  and  detection  techniques  with  the  constraint  to  maximize  use  of  limited 
resources.  Washburn  (2002)  provides  many  mathematical  techniques  to  determine  the  proba¬ 
bility  of  detecting  a  target.  Rozen  (2009)  develops  a  stochastic  dynamic  programming  model  to 
optimize  the  operational  policy  of  an  interdiction  force  performing  a  maritime  interdiction  mis¬ 
sion.  These  works  consider  cases  of  varying  intelligence  about  a  target,  where  best  to  search, 
and  when  to  perform  interdiction,  while  this  thesis  develops  a  grid  of  probabilities  overlaid  on 
the  AOR  and  then  determines  the  target’s  track.  Our  focus  develops  tools  to  generate  the  inputs 
for  search  models  such  as  Pietz  and  Royset  (2013). 

2.2  Piracy  Models 

Hansen  et  al.  (2011)  develop  a  probability  model  to  aid  in  the  location  of  pirates.  Their  model 
dynamically  combines  real-time  intelligence  and  weather  information  to  provide  a  visual  pic¬ 
ture  of  likely  pirate  locations.  Intelligence  such  as  homeports,  vessel  information,  performance 
characteristics,  pirate  activity,  pirate  operations,  and  observations  of  likely  pirate  vessels  are 
inputs  in  the  model.  Additionally,  wind,  wave,  and  current  data  is  supplied  creating  a  "dynam¬ 
ically  coupled"  model  (Hansen  et  al.,  2011).  The  model  produces  a  prediction  of  the  pirate’s 
track  over  time  with  the  constraints  of  intelligence  and  weather.  After  thousands  of  runs  are 
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performed,  the  model  produces  a  probability  map  or  "heat  map"  of  pirate  locations.  This  map 
is  updated  as  additional  intelligence  is  received. 


Building  on  the  Hansen  et  al.  (201 1)  piracy  model.  An  et  al.  (2012)  proposes  an  extension  to  the 
model  that  dynamically  allocates  search  and  interdiction  assets  to  minimize  the  probability  of 
a  successful  pirate  attack.  Their  research  focuses  on  allocating  multiple  search  and  interdiction 
assets  against  multiple  targets  or  pirates.  An  et  al.  (2012)  divide  the  asset  allocation  problem 
into  two  phases.  First,  they  allocate  the  interdiction  assets.  Second,  they  divide  the  area  not 
covered  by  interdiction  assets  into  rectangular  areas  for  search  asset  allocation. 

A  pirate’s  use  of  a  small  vessel  is  significantly  impacted  by  winds  and  sea  state  leading  to  the 
development  of  this  model.  Similarly,  drug  traffickers  operate  small  vessels  during  their  transit 
from  South  America  to  the  major  drug  transit  countries  making  them  susceptible  to  the  same 
types  of  weather.  The  Hansen  et  al.  (2011)  model  entails  a  high  degree  of  complexity  through 
the  simultaneous  analysis  of  weather,  intelligence,  and  pirate  behavior.  This  model  has  been 
extended  to  the  drug  trafficker  problem  but  retains  its  complexity.  Our  goal  is  to  develop  a 
much  simpler  model.  We  analyze  straight  line  paths  between  an  origin  and  destination  without 
the  addition  of  weather  information.  Expanding  this  approach  we  apply  similar  techniques  to 
target  paths  that  involve  multiple  courses  and  non-constant  velocity.  Our  approach  focuses  on 
determining  a  single  target’s  path  through  the  AOR. 


2.3  Probability  Model  -  Optimization  Model  Link 

Pfeiff  (2009),  Gift  (2010),  and  Bessman  (2010)  develop  attacker  -  defender  type  interdiction 
models  to  maximize  the  searcher’s  probability  of  detecting  an  illicit  drug  trafficker.  These 
theses  make  use  of  intelligent  smugglers  who  are  able  to  learn  and  adapt  their  voyage  to  the 
searcher’s  actions.  This  thesis,  while  retaining  some  elements  of  search  and  detection,  is  pri¬ 
marily  concerned  with  determining  the  path  of  the  target  to  include  start  time,  origin,  speed, 
course,  waypoints  (if  applicable),  end  time,  and  destination. 

Pietz  and  Royset  (2013)  and  Pietz  (2013)  develop  an  optimization  model  that  routes  JIATF 
South  detection  and  interdiction  assets  to  maximize  the  amount  of  drugs  interdicted.  This  model 
does  not  take  a  probability  map  directly  as  input,  but  takes  piecewise  linear  paths  to  describe 
the  target  track  and  an  Area  of  Uncertainty  (AOU)  surrounding  the  path.  JIATF  South  planners 
may  not  have  access  to  the  necessary  parameters  to  run  the  optimization  model  and  if  weather 
significantly  impacts  the  target’s  course,  the  inputs  will  be  insufficient.  Our  focus  on  transform- 
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ing  the  probability  map  into  a  path  links  the  probabilities  generated  from  gathered  intelligence 
to  the  optimization  model  of  JIATF  South  assets.  Furthermore,  we  develop  an  algorithm  to 
determine  the  optimal  AOU  surrounding  the  target.  These  efforts  directly  support  the  inputs 
needed  for  the  optimization  model. 

2.4  Statistics  Problem 

The  probability  model  produces  a  set  of  probabilities  indexed  over  time  and  identified  by  the  lat¬ 
itude  and  longitude  of  the  respective  grid  cell  center  point.  To  determine  the  target’s  track  during 
the  time  period,  several  statistical  methods  are  considered.  Principle  curves  are  smooth  curves 
that  pass  through  the  middle  of  the  data  cloud  minimizing  the  distance  from  the  data  (Hastie, 
1984).  These  smooth  curves  give  a  general  idea  of  the  target  path  but  lack  information  needed 
on  specific  waypoints.  Smoothing  splines  provide  similar  descriptions  of  the  target  track  as  in 
principle  curves  but  also  lack  the  needed  granularity  on  when  a  target  changes  course  (Hastie, 
Tibshirani,  &  Friedman,  2009).  Hastie  et  al.  (2009)  discusses  the  complexity  of  determining 
the  appropriate  location  of  knots,  and  suggests  the  use  of  Multivariate  Adaptive  Regression 
Splines  (MARS). 

MARS  employs  a  greedy  algorithm  and  various  approximations  to  reduce  the  complexity  of 
knot  selection  (Hastie  et  al.,  2009).  MARS  has  been  employed  in  areas  of  research  con¬ 
cerning  survival  analysis,  time  series,  ecology,  and  economics  since  its  introduction  by  Fried¬ 
man  (1991).  MARS  is  a  non-parametric  regression  technique  that  automatically  models  non- 
linearities  and  selects  variables  based  on  their  contribution  to  overall  model  fit.  This  results  in 
a  piecewise  linear  model  (Hastie  et  al.,  2009).  For  this  thesis,  we  desire  to  define  a  target’s 
track  as  a  set  of  piecewise  linear  paths  which  is  what  MARS  produces  automatically.  MARS 
allows  paths  with  multiple  courses  or  speeds  to  be  modeled  accurately  by  automatically  de¬ 
termining  locations  of  potential  course  or  speed  changes.  Automatic  selection  of  these  points 
makes  this  technique  preferred  over  other  regression  techniques.  Once  MARS  determines  the 
points  of  course  or  speed  changes,  all  points  are  checked  for  validity  and  any  spurious  points  are 
removed.  MARS  supplies  the  points  of  speed  and  heading  change  and  the  algorithm  builds  an 
ordinary  least  squares  regression  model  utilizing  the  selected  points  of  course  or  speed  change. 

2.5  Area  of  Uncertainty 

When  working  with  probabilities,  there  is  an  inherent  uncertainty  involved.  In  our  case,  this 
uncertainty  resides  in  the  location  of  the  target.  At  each  time  step,  a  set  of  probabilities  are 
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created  over  a  particular  area.  An  AOU  is  calculated  for  each  time  period  using  several  methods, 
which  include  minimum  bounding  box,  convex  hull,  and  a  greedy  algorithm. 

Minimum  bounding  boxes  and  convex  hulls  are  used  in  a  variety  of  applications,  including  col¬ 
lision  avoidance,  hidden  object  determination,  and  shape  analysis.  Johnson  and  Cohen  (1998), 
Lempitsky,  Kohli,  Rother,  and  Sharp  (2009),  and  Huebner,  Ruthotto,  and  Kragic  (2008)  discuss 
these  applications  in  detail.  For  our  purposes,  we  present  a  two-dimensional  set  of  points  that 
we  wish  to  enclose  with  the  minimum  bounding  box.  The  area  of  the  minimum  bounding  box 
then  represents  the  maximum  AOU  of  the  target  at  a  particular  time  and  encloses  the  entire 
probability  mass.  The  convex  hull  further  reduces  the  area  while  still  capturing  all  of  the  prob¬ 
ability  mass.  When  the  maximum  AOU  or  entire  probability  mass  is  not  necessary,  we  propose 
the  use  of  a  greedy  algorithm  to  determine  the  AOU  for  a  user  specified  level  of  probability. 
This  approach  allows  the  user  to  select  the  minimum  amount  of  probability  the  AOU  should 
contain. 
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CHAPTER  3: 
Probability  Model 


JIATF  South  processes  various  amounts  of  intelligence  information  and  makes  decisions  on 
where  and  when  to  employ  available  assets.  Given  this  intel,  we  develop  probability  models 
to  generate  a  heat  map  to  provide  the  operators  with  better  situational  awareness  about  the 
likely  locations  of  targets,  both  in  the  current  period  and  future  periods.  We  first  build  a  simple 
probability  model  that  accounts  for  uncertainty  only  in  the  departure  and  arrival  locations.  From 
this  model,  we  further  enhance  its  functionality  to  account  for  uncertain  speed  and  time  of 
departure  as  well  as  consideration  for  multiple  legs.  We  first  present  the  simplest  model  and 
then  describe  the  embellished  model. 


3.1  Fundamental  Probability  Model 

The  fundamental  one  target  analytic  model  accounts  for  uncertainty  only  in  the  departure  and 
arrival  locations  on  a  rectangular  area  of  operation.  This  setup  does  not  consider  the  geography 
of  realistic  shorelines  and  we  assume  the  target  travels  from  the  bottom  of  the  rectangle  to  the 
top  of  the  rectangle  in  a  straight  line  as  shown  in  Figure  3.1.  We  further  assume  the  target’s 
departure  time  and  velocity  are  known.  We  set  up  the  AOR  by  defining  an  area  from  0  to 
Xmax  units  in  width  and  0  to  Ymax  units  in  height.  Within  the  Xmax  by  Ymax  area,  we  discretize 
the  AOR  into  a  set  of  square  grid  cells  for  future  probability  calculations.  We  further  assume 
that  the  target’s  y  coordinate  at  the  origin  is  0  and  the  y  coordinate  at  the  destination  is  Ymax. 
Figure  3.1  illustrates  the  model  setup. 

The  random  variables  in  this  case  apply  to  the  departure  and  arrival  locations  in  the  x  direc¬ 
tion  and  are  defined  as  D  and  A  respectively.  The  target  velocity  and  departure  time  remain 
constant  and  are  deterministic  in  this  simplified  model.  Additionally,  given  the  departure  and 
arrival  locations,  we  assume  the  target  moves  in  a  straight  line  between  the  locations.  Let  v 
represent  the  target  velocity  and  to  represent  the  target  departure  time.  The  random  variable  D 
follows  a  probability  distribution  Fq  on  the  interval  [D/,  D„]  and  the  random  variable  A  follows 
a  probability  distribution  Fa  on  the  interval  [A/,A„].  Let  a  represent  a  realized  value  of  A  and  d 
represent  a  realized  value  of  D.  We  desire  to  find  the  probability  that  the  target  is  in  a  specific 
box  at  a  specific  time.  We  define  the  box  by  upper  and  lower  coordinates  in  both  the  x  and  y 
directions.  Let  xj  and  xu  define  the  lower  and  upper  edges  of  the  box  in  the  x  direction  and  y/ 
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Simple  Probability  Model  Layout 


x 

Figure  3.1:  Initial  layout  of  the  fundamental  probability  model,  which  shows  the  various  depar¬ 
ture  and  arrival  locations,  target  track,  and  box  of  interest.  The  target  travels  from  the  lower  left 
comer  of  the  area  of  interest  to  the  upper  right  comer  with  constant  velocity. 


and  yu  define  the  lower  and  upper  edges  of  the  box  in  the  y  direction.  Let  (xt,yt)  be  the  coordi¬ 
nates  of  the  target  at  time  t.  The  target  position,  (xt,yt),  is  dependent  on  a,  d,  v,  t,  and  to  and  is 
defined  in  Equations  3.1  and  3.2.  We  establish  an  indicator  variable  equal  to  1  if  the  target  is  in 
the  box  and  0  otherwise.  Equation  3.3  defines  the  indicator  variable.  For  simplicity,  we  define 
square  grid  cells  or  boxes  for  probability  calculations  but  any  arbitrary  shaped  area  can  be  used 
by  modifying  Equation  3.3. 
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(3.1) 


xt(d,a,v,t,tQ ) 


j  |  (a-d)v(t-to) 
V(a-d)2  +  Y^ax 


yt(d,a,v,t,to) 


YmaxV{t  fit) 
\/(a-d)2  +  Y?-ax 


(3.2) 


I{d,a,v,t,to) 


if  xt<x t<  xu  and  yi<yt<  yu 

otherwise 


(3.3) 


Given  a  departure  location,  d  in  [D[,DU]  and  an  arrival  location,  a  in  \Ai,Au\,  we  know  precisely 
whether  the  target  is  in  the  box  at  a  specific  time  from  Equation  3.3.  The  probability  based  on 
an  uncertain  d  and  a  results  from  integrating  over  the  distributions,  Fd  and  Fa,  which  yields 
Equation  3.4. 


r  D;<  rAu 

P[xi  <xt  <xu,yi  <yt  <yu\=  /  I{d,a,v,t,t0)fD(d)fA(a)ddda  (3.4) 

Jd,  JAi 

3.2  Numerical  Illustrations 

After  deriving  a  method  to  determine  the  probability  a  target  is  in  one  particular  box  at  a  specific 
time,  we  expand  this  to  each  time  period  and  all  of  the  boxes  comprising  the  operating  area. 
Specifically,  for  each  time  period,  we  want  to  compute  the  probability  a  target  is  in  any  grid  cell 
making  up  the  operating  area.  We  present  a  numerical  example  to  illustrate  the  application  of 
the  model  over  a  multiple  hour  time  period  covering  the  entire  area  of  operations. 

The  simple  one  target  analytic  model  as  implemented  in  R,  uses  discrete  departure  and  arrival 
distributions.  We  divide  the  departure  location  into  n d  bins  and  the  arrival  location  into  na  bins. 
For  example,  the  uniform  departure  distribution  takes  the  value 

d\  =  D/  +  (Du  —  DA  — — —  with  probability  —  for  1  <  i  <  nci  (3.5) 

na~  1  nd 

and  the  arrival  random  variable,  A,  takes  on  similar  quantities  defined  as  aj.  Thus,  the  integrals 
in  Equation  3.4  become  sums  as  shown  in  Equation  3.6. 
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nd  na 

P[xi  <xt<  xu,yi  <yt<  yu]  =  IE  I(di,aj,v,t,xi,yi,xu,yu)fD(di)fA(aj)  (3.6) 

(=1  j=i 

For  example,  let  D  and  A  be  random  variables  with  uniform  probability  distributions  on  the 
intervals  [Di,Du\  and  [A/,A„]  respectively.  Then  Equation  3.6  becomes 

^  ^ d  fta 

P[x,  <xt  <xUlyi<yt  <yu]  = - £  £ /(d;,ay-,  (3.7) 

nd  na  i=i  j=i 

The  probability  model  produces  target  tracks  from  each  d,  to  each  aj  and  plots  the  associated 
tracks.  When  the  indicator  variables  are  summed  and  multiplied  by  the  associated  density  func¬ 
tions,  we  produce  the  probability  that  the  target  is  in  a  specific  grid  cell  for  each  time  period.  By 
repeating  the  calculations  for  each  time  period  over  the  entire  length  of  the  journey,  we  produce 
a  matrix  of  probabilities.  This  matrix  of  probabilities  is  then  used  to  generate  a  spatial  temporal 
probability  map  for  each  time  period.  An  example  of  the  heat  map  is  shown  in  Figure  3.2  where 
the  departure  location  follows  a  uniform  distribution,  Fq,  on  the  interval  [10,30]  and  the  arrival 
location  follows  a  uniform  distribution,  Fa,  on  the  interval  [70,90]. 


3.3  Enhanced  Probability  Model 

The  addition  of  uncertainty  in  velocity  and  departure  time  further  enhances  the  fundamental 
probability  model.  The  additional  uncertainty  introduces  two  more  random  variables.  We  de¬ 
fine  these  as  V  for  velocity  and  T  for  departure  time.  The  random  variable  V  follows  a  prob¬ 
ability  distribution  Fy  on  the  interval  [ V/ ,  Vu]  and  the  random  variable  T  follows  a  probability 
distribution  Ft  on  the  interval  [7],  Tu\.  Equation  3.4  gains  two  additional  integrals  to  account  for 
the  additional  uncertainty  as  shown  in  Equation  3.8. 


P[xi  <xt<  xu,yi  <yt<  yu ]  = 

rT,<  rV,i  rDu  rAu 

/  /  /  /  I (di,  cij,  v.t,xi,yi,xu, yu)fT{t)fv ( v)fn (d)fA  (a)dtdvddda  (3.8) 

Jt,  Jv,  Jd,  Ja, 

Furthermore,  the  enhanced  probability  model  can  be  used  to  generate  probabilities  for  targets 
with  multiple  legs.  The  legs  are  divided  by  waypoints  or  geographic  positions  where  the  target 
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Figure  3.2:  Heat  maps  from  several  periods  of  the  one  target  simple  analytic  model  utilizing 
uniform  departure  and  arrival  distributions  on  the  intervals  [10,30]  and  [70,90],  respectively. 

changes  course,  speed,  or  both.  This  functionality  of  the  enhanced  model  is  an  extension  of 
the  previous  model  where  rather  than  having  just  two  points  on  the  path,  departure  and  arrival, 
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there  can  can  be  an  arbitrary  number  of  waypoints.  If  there  are  M  total  waypoints  in  addition  to 
the  departure  and  arrival  points,  then  there  will  be  M  additional  integrals. 

As  in  the  fundamental  probability  model,  the  basic  unit  considered  by  the  enhanced  model 
is  probability  and  if  the  mass  of  the  target  is  all  on  the  water,  the  sum  of  the  mass  is  one.  The 
enhanced  model  still  utilizes  square  grid  cells  for  probability  calculations  but  the  area  of  interest 
considers  the  actual  geography  of  the  eastern  Pacific  and  Caribbean.  Figure  3.3  shows  the 
enhanced  model’s  area  of  interest.  By  considering  realistic  geographies  Equation  3.5  becomes 
more  complex  in  that  the  current  1 -dimensional  equations  become  2-dimensional  equations. 


Figure  3.3:  Depiction  of  the  area  of  interest  in  the  enhanced  probability  model.  The  enhanced 
probability  model  considers  the  realistic  geography  of  the  eastern  Pacific  and  Caribbean. 


We  construct  a  model  that  generates  target  location  probabilities  over  the  course  of  the  target’s 
voyage.  The  model  takes  as  inputs:  a  starting  location,  a  destination,  target  speed,  waypoints, 
expected  drug  load,  and  intel  certainty.  The  model  applies  a  predetermined  probability  distribu¬ 
tion  to  the  locations,  speed,  and  departure  time  inputs.  The  user  can  choose  either  a  determinis¬ 
tic  distribution,  uniform  distribution,  or  triangular  distribution.  The  additional  considerations  in 
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the  enhanced  model  introduce  more  complexity  than  the  basic  model  in  Section  3.1  but  the  core 
framework  of  the  model  is  the  same.  Numerically,  implementation  of  the  enhanced  model  only 
requires  additional  for  loops  and  retains  its  simplicity  when  compared  to  the  simulation  models 
such  as  the  piracy  model  discussed  in  Section  2.2 

The  enhanced  model  provides  several  generalizations  to  produce  a  more  realistic  heat  map.  The 
uncertainty  in  velocity  and  departure  time  simply  add  two  integrals  to  Equation  3.4  and  the 
more  complex  area  of  interest  requires  a  more  careful  bookkeeping  of  departure  and  arrival 
locations.  By  including  the  expected  drug  load  as  an  input,  the  model  can  produce  a  heat  map 
in  terms  of  pounds  of  drugs  vice  probability  by  simply  multiplying  the  various  probabilities 
by  the  expected  drug  load.  Considerations  for  target  existence  can  be  made  by  multiplying  the 
probabilities  by  the  intel  certainty.  We  can  model  one  target  or  several  targets  simultaneously  on 
different  tracks  displaying  all  targets  on  one  heat  map,  which  would  allow  JIATF  South  planners 
to  analyze  more  complex  target  scenarios.  Figure  3.4  shows  a  scenario  with  four  targets  on  the 
water  en  route  to  various  destinations  in  Mexico  and  Central  America.  Figure  3.4a  depicts  the 
four  targets  15  hours  into  the  scenario  and  Figure  3.4b  shows  three  targets,  as  one  of  the  four  has 
reached  it’s  destination.  The  introduction  of  departure  time  uncertainty  is  evident  in  the  long 
blue  channel  in  Figure  3.4b.  This  particular  target  has  a  departure  window  of  48  hours  thereby 
spreading  the  probability  mass  across  the  entire  voyage.  We  refer  to  these  cases  as  "channel" 
cases.  The  optimization  model  handles  these  cases  differently  and  requires  additional  inputs. 
We  currently  exclude  these  cases  from  analysis.  Following  the  construction  of  a  more  general 
probability  model,  we  begin  the  process  of  analyzing  the  probability  maps  to  reproduce  the 
target  voyage  characteristics. 
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Figure  3.4:  A  multiple  target  scenario  produced  in  the  enhanced  probability  model.  The  more 
general  probability  model  is  able  to  consider  several  targets  simultaneously,  allowing  JIATF 
South  planners  to  analyze  more  complex  target  scenarios.  The  long  "channel"  is  a  special  case 
when  the  departure  window  is  large.  The  optimization  model  handles  these  differently  and  we 
exclude  these  types  of  cases  from  analysis. 
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CHAPTER  4: 
Path-finding  Model 


A  strong  focus  of  this  research  analyzes  how  to  create  the  target’s  track  from  the  probability  map 
generated  by  the  enhanced  probability  model.  Specifically,  we  determine  the  target’s  earliest 
and  latest  departure  times,  earliest  and  latest  arrival  times,  speed  on  each  leg,  course  on  each  leg, 
location  of  waypoints,  time  of  change  in  course  or  speed,  Area  of  Uncertainty  (AOU),  departure 
location,  and  arrival  location.  This  data  is  then  fed  into  the  optimization  model  described  in 
Section  2.3  to  determine  optimum  employment  of  JIATF  South  assets  to  maximize  the  expected 
amount  of  drugs  seized.  Additionally,  the  analytic  probability  model  provides  the  ability  to 
verify  more  complex  simulation  models  such  as  An  et  al.  (2012)’s  extension  of  piracy  model  to 
counter  drug  operations.  The  path-finding  model  assumes  that  the  probability  data  is  the  only 
data  available  where  p(x,y,t )  is  the  probability  that  a  target  is  located  in  a  grid  cell  centered  at 
(jc,y)  at  time  t.  From  this  data,  we  derive  all  other  parameters  for  the  optimization  model. 

4.1  Optimization  Model 

The  Pietz  and  Royset  (2013)  optimization  model  dictates  the  output  required  from  our  path¬ 
finding  model.  This  optimization  model  is  a  path-constrained  optimal  search  problem  in  con¬ 
tinuous  space  and  time  (Pietz  &  Royset,  2013).  The  model  determines  the  optimal  routing  of 
search  assets  in  the  area  of  operations  to  maximize  the  expected  drug  load.  The  area  of  oper¬ 
ations  is  modeled  as  a  transportation  network,  in  which  the  nodes  are  search  regions  and  the 
arcs  are  the  searcher’s  path  between  moving  search  regions.  These  search  regions  are  areas  of 
uncertainty  concerning  the  target’s  location,  thus  knowledge  of  target  parameters  is  necessary 
to  run  the  optimization  model.  The  inputs  associated  with  the  targets  necessary  to  run  the  opti¬ 
mization  model  are  listed  in  Table  4. 1  and  they  correspond  to  those  defined  in  Pietz  and  Royset 
(2013). 

The  optimization  model  assumes  that  the  targets  travel  in  a  piecewise  linear  motion  at  a  known 
constant  speed  (Pietz  &  Royset,  2013).  The  path-finding  model  derives  the  uncertainty  data 
listed  in  Table  4.1  for  input  into  the  optimization  model.  Initially  we  identify  any  waypoints, 
which  correspond  to  changes  in  a  target’s  course  or  speed  and  divide  the  path  into  the  appro¬ 
priate  number  of  legs.  Subsequent  analysis  determines  velocity,  initial  position,  final  position, 
expected  departure  time,  departure  time  uncertainty,  departure/arrival  position  uncertainty,  and 
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Expected  departure  time  of  target  j  Xj 

Time  uncertainty  range  of  target  j  Xj 

Expected  departure  location  of  target  j  pj 

Expected  arrival  location  of  target  j  pj 

Departure/arrival  location  uncertainty  range  of  target  j  pj 
Expected  value  of  detecting  target  j  qj 

Speed  of  target  j  Uj 

Table  4.1:  The  inputs  required  to  run  the  optimization  model  creating  a  search  vehicle  allocation 

plan.  (From  Pietz  and  Royset,  2013). 


AOET  for  each  leg  of  the  target’s  path.  It  is  possible  some  of  the  uncertainty  data  in  Table  4. 1 
is  directly  available  from  JIATF  South  databases  or  intelligence.  However,  in  a  worst  case  sce¬ 
nario,  planners  have  only  the  probability  map  and  must  determine  the  aforementioned  inputs 
solely  from  this  map  to  produce  a  search  vehicle  allocation  plan  via  the  optimization  model. 

4.2  General  Path-finding  Model  Approach 

The  path-finding  model  approximates  the  path  taken  by  the  target  as  a  sequence  of  legs  where 
the  collection  of  all  legs  produces  the  target’s  entire  track  from  departure  location  to  arrival 
location.  We  define  a  leg  as  a  segment  of  the  target’s  track  that  has  a  constant  heading,  speed, 
and  AOU.  This  creates  a  piecewise  linear  path  for  the  target.  Points  in  time  in  which  the  target’s 
speed,  heading,  or  AOU  changes  signal  the  ending  of  one  leg  and  the  beginning  of  the  next.  We 
address  new  legs  based  on  changes  in  the  AOU  once  all  the  constant  speed-heading  legs  are 
identified.  Figure  4.1  shows  a  potential  path  from  departure  location  to  arrival  location  for  one 
target.  In  this  case,  the  target’s  path  consists  of  two  legs  in  which  the  target  changes  heading 
at  the  waypoint.  The  mechanics  of  the  path-finding  model  automatically  identifies  where  and 
when  a  target  changes  heading  or  speed.  This  allows  us  to  generate  the  piecewise  linear  path 
required  in  the  optimization  model. 

On  each  leg  we  must  produce  linear  tracks  with  a  constant  speed,  heading,  and  AOU.  We  deter¬ 
mine  these  inputs  for  each  target  separately  and  compile  the  inputs  by  leg  for  the  optimization 
model.  We  apply  regression  techniques  to  generate  piecewise  linear  paths  through  the  proba¬ 
bility  mass.  Specifically,  the  path-finding  model  produces  two  linear  regression  equations  of 
the  target’s  location,  one  for  longitude  and  one  for  latitude.  The  coefficients  associated  with  the 
two  linear  regression  equations  are  then  used  to  calculate  the  required  quantities.  The  associ¬ 
ated  mathematical  equations  are  discussed  further  in  Section  4.4  and  numerical  illustrations  are 
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Figure  4.1:  Illustration  of  one  target  traveling  along  two  legs  from  departure  to  arrival  location. 
This  track  consists  of  two  legs  with  a  heading  change  at  the  waypoint. 


provided  in  Sections  4.5  and  4.6. 

From  the  probability  model,  we  analyze  a  set  of  non- zero  probabilities  that  are  spread  over  a 
cellular  grid  in  the  AOR.  For  each  time  period,  the  probabilities  correspond  to  the  center  of  a 
grid  cell  by  their  associated  latitude  and  longitude.  We  define  n  as  the  number  of  grid  cells  in 
the  AOR,  pi(t)  as  the  probability  the  target  is  in  cell  i  at  time  t,  and  (jc/.y,-)  as  the  center  point 
of  grid  cell  i.  Figure  3.2  illustrates  four  different  hours  of  the  probability  mass  in  the  AOR  with 
each  grid  cell  colored  based  on  the  target’s  probability  of  being  in  that  cell  at  the  associated 
time.  This  probability  mass  changes  size  and  shape  and  moves  over  time  as  the  target  travels 
from  the  departure  location  to  the  arrival  location.  At  each  time  step,  we  calculate  the  expected 
value  of  the  target’s  longitude  or  x  component  where  E  [X]  —  Y!!=i  Pi(t)xi  and  the  expected  value 
of  the  target’s  latitude  or  y  component  where  E\Y]  =  £”=1  p, 

Since  the  probability  mass  is  spread  over  several  grid  cells  during  each  time  step,  weighted 
target  tracks  could  be  produced  but  for  simplicity  and  limitations  of  the  techniques  employed, 
we  use  expected  values. 

The  resultant  expected  x  and  y  coordinates  are  then  plotted  in  Figure  4.2  for  a  visual  picture  of 
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the  target’s  expected  track.  Figure  4.2  illustrates  the  same  track  from  Figure  4.1  but  uses  the 
expected  values  of  X  and  Y.  As  in  Figure  4.1,  Figure  4.2  shows  two  legs  with  a  heading  change 
near  86°W  3°N  (274°,  3°). 


Expected  Target  Track 


Figure  4.2:  Plot  of  one  target’s  expected  track  with  two  legs  in  the  AOR. 

Based  on  the  plot  of  latitude  versus  longitude,  we  can  easily  see  a  change  in  heading  but  dis¬ 
cerning  departure  and  arrival  locations  or  speed  changes  is  not  possible.  The  necessary  time 
element  is  not  depicted  in  Figure  4.2.  This  leads  us  to  plot  longitude  versus  time  and  latitude 
versus  time  as  shown  in  Figures  4.3  and  4.4  respectively.  These  plots  become  the  primary  means 
of  determining  the  necessary  target  parameters  in  Table  4.1. 

Our  regression  analysis  utilizes  the  probability  data  between  the  latest  departure  time  and  the 
earliest  arrival  time,  which  corresponds  to  the  time  that  the  entire  probability  mass  is  on  the 
water.  These  points  are  specifically  highlighted  in  Figures  4.3  and  4.4.  Any  change  in  slope 
between  these  points  signifies  a  change  in  the  target’s  speed,  heading,  or  possibly  both.  Changes 
in  slope  outside  this  interval  may  produce  spurious  changes  in  course  or  speed  due  to  endpoint 
effects  and  are  not  important  for  optimization.  Figures  4.3  and  4.4  both  show  a  change  in  slope 
around  hour  45,  signifying  a  change  in  the  target’s  motion.  The  points  where  a  change  in  slope 
exists  in  either  the  x  direction  or  y  direction  indicates  where  one  leg  ends  and  the  next  leg 
begins. 

Detecting  two  legs  is  obvious  in  Figure  4.1  because  there  is  a  considerable  change  in  course. 
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Expected  Target  Track  in  Degrees  Longitude  Over  Time 


Figure  4.3:  Plot  of  one  target’s  expected  track  with  two  legs  in  degrees  longitude  over  time. 


Expected  Target  Track  in  Degrees  Latitude  Over  Time 


Figure  4.4:  Plot  of  one  target’s  expected  track  with  two  legs  in  degrees  latitude  over  time. 


However,  a  change  in  speed  with  a  constant  course  can  not  be  determined  with  a  static  figure. 
Consider  the  example  in  Figure  4.5  where  the  target  is  traveling  in  a  straight  line  from  the 
departure  location  to  the  arrival  location. 

Based  on  Figure  4.5,  it  appears  one  leg  exists  in  this  target’s  path.  Looking  at  the  expected  track 
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Figure  4.5:  Illustration  of  a  potential  target  traveling  in  a  straight  line  path  from  the  departure 
location  to  the  arrival  location. 

in  Figure  4.6  we  are  inclined  to  the  same  conclusion. 

Expected  Target  Track 
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Figure  4.6:  Plot  of  the  one  target  straight  line  path  in  degrees  latitude  versus  degrees  longitude. 

A  closer  look  at  this  figure  shows  the  spacing  between  expected  locations  increasing.  This 
increase  in  space  between  expected  locations  signifies  the  target  has  changed  speed.  We  are  still 
missing  the  time  element  and  cannot  tell  whether  the  target  is  traveling  northwest  or  southeast 
but  we  can  tell  speed  has  either  increased  or  decreased  during  the  journey.  When  we  consider 
the  plots  in  Figures  4.7  and  4.8,  it  is  evident  the  target  is  traveling  northwest  and  increases  speed 
during  the  journey  based  on  the  increase  in  slope  at  hour  25. 
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Expected  Target  Track  in  Degrees  Longitude  Over  Time 
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Figure  4.7:  Plot  of  the  one  target  straight  line  path  in  degrees  longitude  over  time. 


Expected  Target  Track  in  Degrees  Latitude  Over  Time 
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Figure  4.8:  Plot  of  the  one  target  straight  line  path  in  degrees  latitude  over  time. 


This  case  clearly  illustrates  the  necessity  of  analyzing  the  x  and  y  components  of  the  target’s 
motion  separately  and  the  need  to  employ  a  technique  capable  of  producing  piecewise  linear 
equations. 


31 


4.3  Path-finding  Model 


We  desire  to  provide  JIATF  South  with  a  path-finding  model  capable  of  performing  this  analysis 
automatically.  To  meet  this  need  we  employ  a  technique  known  as  MARS.  MARS,  as  described 
in  Section  2.4  automatically  fits  a  piecewise  linear  regression  model  to  the  probability  data  by 
automatically  selecting  "knots"  or  points  of  change  in  course  or  speed. 

MARS  uses  a  forward  model  building  and  backward  term  deletion  procedure  to  produce  the 
final  model  (Hastie  et  al.,  2009).  During  the  forward  model  building  procedure,  MARS  uses  a 
set  of  basis  functions  as  shown  in  Equations  4.1  and  4.2  as  input. 


0  -t)+  = 


if  x>t, 
otherwise 


(4.1) 


(t~x)+  = 


if  x<t, 
otherwise. 


(4.2) 


For  each  input,  MARS  forms  a  set  of  reflected  pairs,  (Xj  — 1)+  and  (/  —  Xj)  . ,  for  each  value  of 
that  input  (Hastie  et  al.,  2009).  The  resulting  model  is  of  the  form 


M  N 

f(X)  =  Po+EPm(X-tm)++EPn(tn-X)+  (4.3) 

m=  1  ;z=l 

where  tm  and  tn  are  the  candidate  times  of  course  or  speed  changes.  Additionally,  the  model 
form  can  include  a  product  of  two  or  more  basis  functions  (Hastie  et  al.,  2009).  For  a  particular 
basis  function,  MARS  estimates  the  associated  coefficient,  /3,  by  minimizing  the  residual  sum- 
of- squares  (Hastie  et  al.,  2009).  The  core  functionality  of  MARS  lies  in  the  selection  of  the 
basis  functions.  Each  step  through  the  forward  model  building  process  a  new  basis  function 
pair  and  its  products  of  all  functions  in  the  model  are  evaluated  for  inclusion  in  the  model.  The 
term  that  produces  the  largest  decrease  in  error  is  then  added  to  the  model.  For  a  full  description 
of  MARS  see  Friedman  (1991)  and  Hastie  et  al.  (2009). 

Consider  the  target  scenarios  presented  in  Figures  4.1  and  4.5.  As  previously  discussed,  these 
targets  follow  tracks  that  consists  of  two  legs  identified  by  a  course  change  in  the  expected  track 
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and  a  speed  change  respectively.  Utilizing  MARS,  the  path-finding  model  automatically  locates 
changes  in  slope  as  seen  in  Figures  4.9  through  4.1 1  and  produces  final  fitted  models  of  the  form 


x  —  /3o  +  /3i(t  —  45)+  +  ^2(45  —  f)+ 


(4.4) 


y  =  j80  +  ft  (f  -  37)+  +  k(t  -  49)+  +  /33(49  -  f)+.  (4.5) 


Expected  Track  in  Degrees  Longitude 
Fith  with  MARS  and  Linear  Regression 


Figure  4.9:  Plot  of  one  target’s  track  with  two  legs  in  degrees  longitude  over  time.  The  resulting 
MARS  model  fits  the  path  extremely  well,  whereas  ordinary  least  squares  regression  is  unable 
to  automatically  identify  the  course  change. 

Figures  4.9  through  4.11  show  the  resulting  MARS  model  fit  and  illustrate  the  poor  resulting 
fit  by  using  standard  linear  regression.  Standard  regression  does  not  adequately  identify  the 
slope  changes,  thus  it  does  not  produce  the  piecewise  linear  track  that  MARS  does.  Note, 
Figure  4.10  displays  two  slope  changes  and  Equation  4.5  includes  different  basis  functions  than 
Equation  4.4.  In  fitting  models  in  the  x  and  y  directions  separately,  MARS  produces  a  best  fit 
model  irregardless  of  the  number  of  basis  functions  in  the  other  direction.  A  few  reasons  exist 
for  the  additional  knot  and  are  discussed  further  in  Section  4.6  as  well  as  the  case  parameters 
and  associated  model  output. 
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Expected  Track  in  Degrees  Latitude 
Fit  with  MARS  and  Linear  Regression 


Figure  4.10:  Plot  of  one  target’s  track  with  two  legs  in  degrees  latitude  over  time.  The  resulting 
MARS  model  fits  the  path  extremely  well,  whereas  ordinary  least  squares  regression  is  unable 
to  automatically  identify  the  course  change. 


Expected  Target  Path 
Fit  with  MARS  and  Linear  Regression 


Figure  4.11:  Plot  of  one  target’s  track  with  two  legs  in  degrees  latitude  versus  degrees  longi¬ 
tude.  The  resulting  MARS  model  fits  the  path  extremely  well,  whereas  ordinary  least  squares 
regression  is  unable  to  automatically  identify  the  course  change. 
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At  this  point,  the  basis  functions  represent  potential  waypoints  along  the  target’s  route.  Next, 
the  path-finding  model  checks  the  validity  of  these  waypoints  by  first  evaluating  the  change 
in  slope  of  reflected  pairs,  and  then  evaluating  the  time  between  waypoints  as  compared  to 
the  time  required  for  the  probability  mass  to  complete  the  associated  turn.  These  concepts  are 
demonstrated  in  greater  detail  in  Sections  4.5  and  4.6 

Once  we  know  which  waypoints  are  valid,  we  construct  a  new  data  set  to  account  for  the  way- 
points.  We  add  a  new  data  field,  wpi ,  for  each  valid  waypoint  where  wpt  equals  the  value 
of  the  associated  basis  function.  If  twp  is  the  time  of  the  respective  waypoint,  then  we  define 
wp  =  (t-twp)+. 

The  path-finding  model  uses  this  data  set  and  ordinary  least  squares  regression  to  produce  the 
final  model,  which  is  of  the  form 


x  — 


i 

Po+Plt  +  Y,PiwPi 


i=  1 


(4.6) 


and 


y  =  Po  +  Pit  +  Y,PiwPi  (4-7) 

i=  1 

where  j  represents  the  number  of  waypoints.  The  additional  data  fields,  identified  by  wpt,  enable 
ordinary  least  squares  regression  to  produce  a  piecewise  linear  fit. 


4.4  Input  Generation 

The  probability  data  and  the  path-finding  model  produce  the  uncertainty  parameters  listed  in 
Table  4.1  to  generate  a  search  vehicle  allocation  plan.  Recall  we  have  uncertainty  in  departure 
location,  arrival  location,  departure  time,  and  velocity  and  in  a  worst  case  scenario  JIATF  South 
planners  will  only  have  the  probability  map  to  generate  the  needed  parameters  for  optimal 
allocation  of  their  search  assets.  The  following  section  describes  how  each  of  the  inputs  are 
generated  from  either  the  probabilities  or  the  model  output.  The  desired  end-state  equates  to  a 
list  of  parameters  for  each  target  and  leg  of  the  associated  target’s  path. 
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4.4.1  Times 

The  introduction  of  T,  the  target’s  random  departure  time,  in  the  enhanced  probability  model 
produces  uncertainties  in  the  departure  time  and  subsequently  the  arrival  time  for  a  target  j.  The 
uncertainty  results  in  a  range  of  possible  departure  times  between  [Ti,Tu]  and  a  corresponding 
range  of  arrival  times.  The  specific  times  for  each  leg  required  in  the  Pietz  and  Royset  (2013) 
optimization  model  are  listed  in  Table  4.2.  Additionally,  we  define  the  earliest  departure  time 
of  target  j  from  the  origin  as  zju  liestdep  and  the  latest  arrival  time  at  the  destination  as  xl-aestarr 
for  computational  purposes. 


Expected  departure  time  of  target  j  Zj 
Time  uncertainty  range  of  target  j  z J 

Latest  departure  time  of  target  j  z"un 

Earliest  arrival  time  of  target  j  z ™ax 

Table  4.2:  Times  required  on  each  leg  for  the  optimization  model  to  generate  a  search  vehicle 
allocation  plan.  (From  Pietz  and  Royset,  2013) 


Pietz  and  Royset  (2013)  discuss  calculating  x'?in  and  z'"ax  from  the  uncertainty  data,  r.  and  f,. 
Our  path-finding  model  first  determines  z e.arhestdeP ,  x"un,  Zj,ax,  and  z l.atestarr  for  the  first  and 
last  leg  respectively  and  then  we  calculate  Zj  and  fj.  It  is  important  to  note  that  these  initial 
time  determinations  are  for  the  first  leg  and  the  last  leg  and  are  determined  directly  from  the 
probability  data.  Each  leg  of  the  target’s  path  requires  the  identification  of  the  times  listed  in 
Table  4.2.  MARS  generates  these  times  for  the  intermediate  waypoints,  whereas  the  probability 
data  generates  the  times  associated  with  the  origin  and  destination.  Figure  4.12  illustrates  the 
associated  times.  In  this  case  more  than  one  leg  exists,  thus  the  intermediate  departure  and 
arrival  times  are  determined  through  regression  analysis. 

The  earliest  departure  time  is  the  first  time  period  that  any  non-zero  probability  exists  whereas 
the  latest  departure  time  is  the  first  time  period  that  the  sum  of  the  probabilities  is  equal  to  1. 
These  are  expressed  in  Equations  4.8  and  4.9  respectively. 


earliestdep 

Zj  =  min  • 


l£p*0)  >° 

i=  1 


Tmin 

Tj 


=  min 


=  1  | 


(4.8) 

(4.9) 
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Expected  Target  Track  in  Degrees  Longitude  Over  Time 


Figure  4.12:  The  initial  time  determinations  in  the  path-finding  model  represent  the  earliest  and 
latest  departure  times  on  the  first  leg  and  the  earliest  and  latest  arrival  times  on  the  last  leg. 


Similarly  for  the  arrival  times,  the  earliest  arrival  time  is  the  maximum  time  that  the  probabilities 
sum  to  1  and  the  latest  arrival  time  is  the  maximum  time  the  probabilities  sum  to  greater  than 
zero.  These  are  expressed  in  Equations  4.10  and  4.1 1  respectively. 


tT  =  maxjt|J>(0 


(4.10) 


-late  Starr 

Lj 


=  max 


AY,Pi{t) 

i=l 


(4.11) 


It  is  helpful  to  view  the  origin  as  a  source  and  the  destination  as  a  sink.  The  probability  mass 
begins  to  grow  from  the  source  and  eventually  100%  of  the  probability  mass  has  left  the  source, 
thus  representing  the  earliest  and  latest  departure  times  respectively.  As  the  target  travels  across 
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the  AOR,  the  probability  mass  begins  to  disappear  into  the  sink  until  all  of  the  mass  has  finally 
collected  in  the  sink,  thus  representing  the  earliest  and  latest  arrival  times.  When  only  one  leg 
exists  for  a  target,  all  times  are  determined  directly  from  the  probability  data.  The  departure 
time  uncertainty  range  of  target  j,  fj,  is  simply  the  difference  between  the  earliest  departure 
time  and  latest  departure  time.  Based  on  xj"hest  ep  and  z"un,  we  determine  Zj  and  Zj  using 
Equations  4.12  and  4.13  respectively. 

^  ^ min ^earliest dep  ^  12) 

J  J  J 


i-  =  —  ^ear^est^eP^j  (4  i3) 

If  the  departure  window  is  large  enough,  the  probability  mass  can  grow  from  the  source  to  the 
sink  without  ever  fully  leaving  the  source.  These  cases  are  referred  to  as  "channel"  cases  and 
the  optimization  model  handles  them  slightly  different.  An  example  is  shown  in  Figure  3.4b. 
At  this  point,  we  only  focus  on  cases  where  the  probability  mass  develops  entirely  on  the  water 
before  beginning  to  absorb  into  the  sink.  For  purposes  of  the  path-finding  model,  we  are  mainly 
concerned  with  the  latest  departure  time  and  the  earliest  arrival  time.  In  other  words,  the  path¬ 
finding  model  only  uses  the  probability  data  when  the  probabilities  sum  to  1,  which  is  the 
probability  data  between  the  first  z1™"  and  the  last  z'jmx. 

4.4.2  Position 

Each  position  is  composed  of  a  horizontal  component  and  a  vertical  component  and  the  random 
variable  D  introduces  uncertainty  in  the  departure  location  on  the  interval  [D{,DU\.  The  regres¬ 
sion  models,  when  run  between  the  first  z'j1"1  and  the  last  z"mx ,  result  in  a  set  of  fitted  values 
between  these  times.  Essentially,  they  are  an  estimate  of  the  target’s  horizontal  and  vertical 
position  in  the  AOR  for  each  time  step  between  departure  and  arrival.  The  optimization  model 
requires  position  inputs  for  each  leg,  specifically  a  starting  position  and  an  ending  position. 
Pietz  and  Royset  (2013)  define  the  origin  of  the  AOR  as  0°  N  (0°)  and  100°  W  (260°)  where  1° 
in  any  direction  equals  60  nm.  All  positions  are  measured  from  this  origin. 

The  fitted  values  from  the  regression  model  directly  correspond  to  the  position  of  target  j. 
Knowing  T™”  and  z'jmx  for  each  leg  allows  us  to  determine  the  target’s  position  by  substi¬ 
tuting  zfm  and  z"tax  into  the  final  models.  The  resulting  (x  z min  5  yzmin  )  and  (. xxmax,yijmax )  are 

J  J  t7  ij  j 
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then  converted  to  nautical  miles  from  the  AOR  origin  where  Pj(x)  —  (xt  —  .ro)  *  60  and  Pjiy)  — 

iyt-yo)  *60. 

Pietz  and  Royset  (2013)  define  Pj  and  p7  as  the  expected  departure  location  and  the  expected 
arrival  location  of  target  j.  In  other  words,  Pj  and  pj  are  the  target’s  position  at  % j  and  the  last 
z'jUlx  respectively.  An  accurate  estimate  of  the  expected  departure  location  is  required  for  the 
optimization  model,  but  an  accurate  estimate  of  the  expected  arrival  location  is  not  required 
as  long  as  its  along  the  same  track.  Since  we  regress  on  the  probability  data  between  the  first 
t'""!  and  the  last  zfax,  our  estimate,  based  solely  on  regression,  of  p,-  is  always  further  from  the 
true  origin.  This  fact  does  not  impact  the  path-finding  model,  however;  to  correct  this,  we  use 
the  vertical  and  horizontal  speeds  to  move  the  position  at  the  first  T'jun  the  distance  of  one-half 
the  departure  time  uncertainty  times  the  respective  speed.  For  example,  with  a  departure  time 
uncertainty  of  4  hours  and  a  vertical  speed  of  10  knots  and  a  horizontal  speed  of  12  knots,  the 
Pj  for  the  first  leg  moves  20  nm  in  the  vertical  direction  and  24  nm  in  the  horizontal  direction 
towards  the  true  origin. 


4.4.3  Speed 

The  optimization  model  assumes  the  searcher  knows  the  speed  of  target  j  on  each  leg  of  the 
path  and  defines  this  speed  as  Uj  (Pietz  &  Royset,  2013).  This  assumption  then  applies  to 
the  entire  leg  without  considering  the  possibility  of  speed  changes  during  the  leg.  Recall,  a 
change  in  speed  or  heading  is  cause  for  the  introduction  of  a  new  leg  since  the  optimization 
assumes  a  constant  speed  on  each  leg.  Considering  the  worst  case  scenario,  where  only  the 
probability  data  exists,  we  desire  to  determine  the  target’s  speed  on  each  leg.  The  coefficients 
of  the  regression  model  in  Equations  4.6  and  4.7  are  used  to  calculate  speed.  The  coefficients 
in  the  linear  models,  excluding  the  intercept  (/3o),  represent  the  speed  in  the  x  and  y  directions. 
Overall  target  speed  is  then  defined  as  U  =  .  / /3A2  +  /3.2. 


In  a  multiple  leg  case,  the  regression  techniques  produce  a  piecewise  linear  model  for  the  entire 
journey  where,  depending  on  the  time,  certain  coefficients  are  active.  Recall  the  MARS  dis¬ 
cussion  in  Section  4.3,  the  basis  functions  equal  zero  when  they  evaluate  to  a  negative  quantity. 
Therefore,  for  times  when  a  basis  function  equals  zero  the  respective  coefficient  does  not  con¬ 
tribute  to  the  target’s  speed.  The  coefficients  only  contribute  to  the  speed  calculation  when  the 
basis  function  evaluates  to  a  non-negative  number.  Knowing  the  T™  and  T"iax  for  each  leg,  the 
path-finding  model  determines  which  coefficients  are  active  and  based  on  the  active  coefficients 
computes  the  speed  for  the  respective  leg. 
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4.4.4  Area  of  Uncertainty 


An  AOU  results  from  uncertainty  in  both  the  departure  time  and  the  departure  location.  This 
area  is  essentially  the  probability  mass  moving  across  the  water.  Just  as  a  constant  speed  on  each 
leg  is  required,  so  is  a  constant  area  along  each  leg.  Changes  in  the  AOU  also  would  indicate 
the  need  for  additional  legs.  Determining  the  area  of  the  probability  mass  as  it  moves  across 
the  water  is  rather  straightforward.  We  can  simply  count  the  number  of  grid  cells  that  contain 
probabilities  greater  than  0.  Knowing  the  size  of  each  grid  cell  and  the  number  of  grid  cells 
with  a  non-zero  probability,  we  simply  multiply  the  two  and  determine  the  associated  AOU. 
This  method  while  straight  forward  does  not  result  in  a  constant  area  over  time.  In  general,  the 
area  changes  every  time  step  and  the  user  may  not  require  100%  coverage  of  the  AOU.  We 
explore  several  methods  to  determine  the  AOU. 

The  first  approach  determines  the  area  of  the  minimum  bounding  box  at  each  time  step  of 
the  target’s  path.  This  method  captures  the  entire  probability  mass  and  results  in  the  largest 
AOU.  Additionally,  it  is  possible  for  the  minimum  bounding  box  to  capture  locations  with  zero 
probability,  as  evident  in  Figure  4.13.  Therefore,  the  second  approach  determines  the  convex 
hull  of  the  probability  mass.  Again,  the  convex  hull  captures  all  of  the  probability  mass  but  with 
minimal  area.  Figure  4.13  depicts  the  minimum  bounding  box  and  convex  hull  at  several  time 
steps.  These  methods  successfully  bound  the  entire  probability  mass  and  capture  100%  of  the 
area. 

The  third  approach  uses  a  greedy  algorithm  to  capture  a  user  defined  amount  of  the  probability. 
Throughout  this  thesis  we  assume  a  probability  threshold  of  90%.  For  each  time  between  the 
first  r jun  and  the  last  z'’iax,  the  greedy  algorithm  selects  the  node  or  cell  with  the  highest  prob¬ 
ability  and  then  creates  a  set  of  cells  that  neighbor  the  initial  cell.  The  algorithm  proceeds  by 
examining  the  set  of  neighbors  and  selects  the  neighbor  with  the  highest  probability  and  adds 
it  to  the  AOU  set.  Next  it  calculates  the  cumulative  probability  of  the  AOU  set  and  proceeds  to 
look  at  the  set  of  neighbors  including  any  new  neighbors,  again  adding  the  neighbor  with  the 
highest  probability.  This  procedure  continues  until  the  cumulative  probability  is  greater  than  or 
equal  to  the  selected  threshold.  Once  we  reach  the  specified  threshold,  we  simply  determine  the 
amount  of  cells  in  the  set  and  multiply  by  the  area  of  a  single  cell  to  calculate  the  total  AOU. 
Figure  4.13  shows  the  grid  cells  that  are  included  in  the  greedy  algorithm’s  AOU  set. 

Figure  4. 14  shows  the  results  of  the  three  methods  over  time.  The  AOU  found  using  the  mini¬ 
mum  bounding  box  tends  to  result  in  a  considerably  higher  area  since  it  encompasses  the  entire 


40 
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Figure  4.13:  Examples  of  the  minimum  bounding  box,  convex  hull,  and  greedy  algorithm 
applied  to  a  target’s  probability  mass  over  several  time  steps  of  a  target’s  path.  The  minimum 
bounding  box  and  convex  hull  encompass  the  entire  probability  mass  at  each  time  step,  while 
the  greedy  algorithm  results  in  the  smallest  AOU. 


Minimum  Bounding  Box 
Convex  Hull 
Greedy  Algorithm 
Max  Probability 
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probability  mass  and  additional  area  with  zero  probability.  The  convex  hull  performs  slightly 
better  but  it  is  still  greater  than  the  greedy  approach.  By  taking  a  probability  threshold  of  90%, 
we  are  able  to  significantly  reduce  the  size  of  the  AOU  while  still  ensuring  we  capture  a  signif¬ 
icant  amount  of  the  probability  mass.  Additionally,  the  greedy  algorithm  allows  for  flexibility 
and  if  desired,  the  user  may  select  a  threshold  of  100%. 


Target  Area  of  Uncertainty 


Figure  4.14:  Comparison  of  the  AOUs  resulting  from  the  three  methods  explored.  The  greedy 
algorithm  consistently  produces  areas  less  than  the  other  two  approaches,  while  the  minimum 
bounding  box  results  in  significantly  higher  areas. 


With  the  AOU  specified  for  each  time  step  of  the  target’s  path,  we  now  desire  to  create  a  piece- 
wise  constant  area  for  the  entire  journey.  The  area  analysis  utilizes  the  functionality  of  MARS 
and  automatically  fits  a  piecewise  linear  model  to  the  area  data.  The  model  smooths  the  area 
data  but  does  not  produce  the  piecewise  constant  areas  required.  Figure  4.15  shows  the  MARS 
model  fit  to  the  areas  from  the  greedy  algorithm.  This  implementation  of  MARS  forces  a  mini¬ 
mum  span  between  knots  equal  to  the  departure  time  uncertainty.  The  minimum  span  prevents 
MARS  from  introducing  knots  at  points  just  before  or  after  a  new  leg. 

Additionally,  for  each  constant  speed  and  heading  leg,  we  build  a  separate  model  since  we 
desire  to  create  the  minimum  number  of  new  legs  based  on  area,  and  we  know  the  number  of 
legs  based  solely  on  speed  and  heading  changes.  Figure  4.15  depicts  the  MARS  models  fit  to 
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MARS  Model  Fit  to  Greedy  AOU 


Time  (Hours) 

Figure  4.15:  A  MARS  model  smooths  the  AOU  over  time.  A  different  model  is  built  for  each 
leg  with  a  constant  speed  and  heading.  The  fitted  values  from  the  model  are  further  analyzed  to 
produce  a  piecewise  constant  AOU. 


each  leg  of  the  target’s  path.  To  determine  the  number  of  additional  legs  based  on  a  change  in  the 
AOU,  we  calculate  the  maximum  number  of  legs  that  we  can  fit  in  each  constant  speed-heading 
leg  based  on  two  different  criteria.  First,  we  use  the  minimum  span  to  determine  the  number  of 
additional  legs  that  can  be  fit  into  each  constant  speed-heading  leg.  This  is  simply  the  length 
of  the  leg  in  time  divide  by  the  departure  time  uncertainty.  The  second  method  takes  a  user 
selected  percent  change  in  area,  which  is  currently  set  at  25%,  and  determines  the  number  of 
times  a  25%  change  in  area  occurs.  Then  we  take  the  minimum  between  the  two  methods.  This 
gives  us  the  number  of  new  legs  to  insert  based  on  area  changes.  We  perform  this  procedure 
on  each  constant  speed-heading  leg.  Considering  the  slope  of  the  MARS  model,  we  assign 
areas  to  the  new  legs  at  the  maximum  area  of  the  new  interval.  Figure  4.16  shows  the  resulting 
piecewise  constant  area. 

In  this  particular  case,  there  are  two  constant  speed-heading  legs  where  the  course  changes  at 
hour  30.  Based  on  the  AOU  analysis,  we  divide  the  first  constant  speed-heading  leg  into  two 
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Piecewise  Constant  AOU 
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Figure  4.16:  Illustration  of  a  piecewise  constant  AOU. 


segments.  The  first  segment  ends  at  hour  21,  at  which  point,  the  AOU  increases.  The  second 
constant  speed-heading  leg  remains  unchanged  because  the  AOU  does  not  change  by  more  than 
25%. 


4.5  Straight  Line  Model  and  Illustration 

Assume  there  is  one  target  of  interest  that  travels  from  the  origin  to  the  destination  with  a 
constant  speed  and  heading.  Figure  4.17  shows  the  target’s  track  from  departure  location  to 
arrival  location.  We  illustrate  the  path-finding  model  by  first  entering  the  case  parameters  into 
the  enhanced  probability  model.  The  case  parameters  as  shown  in  Table  4.3  include  uncertainty 
in  departure  and  arrival  locations  and  departure  time. 

Running  the  probability  model  produces  the  probability  data  over  time  necessary  for  the  path¬ 
finding  model.  Figure  4.18  depicts  several  hours  of  the  target’s  location  as  it  moves  through  the 
AOR.  In  a  worst  case,  this  is  all  the  JIATF  South  planners  have  and  they  desire  to  determine  the 
information  in  Table  4.3. 
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Figure  4.17:  Illustration  of  a  target’s  track  from  departure  to  arrival  location.  The  target  follows 
a  direct  path  at  a  constant  speed  to  illustrate  the  simplest  path-finding  case. 


Approximate  payload 

Speed 

Origin 

Departure  uncertainty 
Destination 
Arrival  uncertainty 
Departure  time 
Departure  time  uncertainty 
Intel  certainty 


5,000  kg 
15  knots 
76.3°  W  9°  N 
±30  nm,  Uniform 
83.5°  W  14°  N 
±30  nm.  Uniform 
33  hours 

±4  hours,  Uniform 
95% 


Table  4.3:  Parameters  for  a  simple  straight  path  from  origin  to  destination  with  constant  velocity. 
Uncertainty  data  is  listed  for  departure  and  arrival  locations  and  departure  time. 


In  this  case,  the  probability  model  generates  the  heat  map  for  78  hours  using  the  associated 
parameters  from  Table  4.3.  In  the  context  of  this  thesis,  hour  1  represents  the  first  simulated 
time  step  and  hour  78  represents  the  last  time  step.  This  results  in  78  probability  data  files,  one 
for  each  time  step,  that  we  read  into  the  path-finding  model.  Each  file  contains  a  probability 
for  each  grid  cell  in  the  AOR  corresponding  to  the  likelihood  the  target  is  in  that  grid  cell  at  the 
associated  time.  The  path-finding  model  analyzes  the  probabilities  to  establish  the  times  listed 
in  Table  4.2  using  Equations  4.8  through  4.13.  Table  4.4  lists  the  resulting  times. 
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(a)  Hour  37 


(b)  Hour  52 


(c)  Hour  64 


Figure  4.18:  A  simple  one  target  scenario  produced  in  the  enhanced  probability  model.  The 
snapshots  are  at  the  latest  departure  time  (hour  37),  an  intermediate  point  (hour  52),  and  at  the 
earliest  arrival  time  (hour  64). 


Once  the  first  departure  time  (t "“")  and  last  arrival  time  (t”mx)  are  known,  the  path-finding 
model  proceeds  by  applying  MARS  to  the  probability  data  between  these  times.  The  resulting 
form  of  the  MARS  model  is 


jcf  =  280.42  -0.205(1  -49)+  +  0.205(49 -t)+  (4.14) 


and 
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Time 

Actual 

Model  Result 

^earliestdep 

29  hours 

29  hours 

First  Tmln 

37  hours 

37  hours 

Fast  xmax 

64  hours 

64  hours 

j late starr 

NA 

72  hours 

f 

8  hours 

8  hours 

T 

33  hours 

33  hours 

Table  4.4:  Case  time  information  from  path-finding  model  as  compared  to  the  actual  inputs  for 
a  target  traveling  along  a  straight  path. 


y,  =  12.13  +  0.142(4  —  55)+  —0.142(55  —  f)+.  (4.15) 

From  Equations  4.14  and  4.15,  we  see  MARS  introduces  the  basis  functions  (t  —  49)+  and  (49  — 
?)+  in  the  x  direction  and  (/  —  55)+  and  (55  —  /)+  in  the  y  direction.  Further  examination  reveals 
the  coefficients  associated  with  the  reflected  pairs  have  equal  magnitude  but  opposite  sign.  This 
is  a  result  of  the  MARS  model  fitting  procedure.  Continuing  the  path-finding  algorithm  with 
this  model  results  in  a  track  composed  of  three  legs  but  the  same  speed  and  heading  on  each  leg. 
In  this  particular  case,  it  is  easy  to  determine  that  these  points  are  not  truly  changes  in  speed  or 
heading.  In  a  more  general  case,  it  is  necessary  to  determine  if  a  point  introduced  by  MARS  is 
really  a  change  in  speed  or  heading. 

We  apply  a  filter  to  the  knots  or  cuts  introduced  by  MARS.  The  filter  first  examines  the  change 
in  slope  at  each  reflected  pair.  If  the  absolute  difference  in  slope  for  both  models  changes  by 
less  than  0.01  the  filter  does  not  add  them  to  the  expected  track  and  path-finding  continues  with 
the  expected  target  locations  as  defined  in  Section  4.2.  In  cases  where  the  reflected  pairs  or 
cuts  are  legitimate,  we  add  an  additional  data  field  to  the  expected  track  data.  This  data  field 
is  identified  by  the  associated  time  and  populated  with  the  result  of  the  basis  function  for  each 
time  step  in  the  path.  For  the  moment  let  us  consider  the  reflected  pair  at  hour  49  as  a  legitimate 
time  to  split  the  journey  in  the  x  direction  into  two  legs.  The  expected  track  data  would  then 
get  an  additional  data  field  containing  the  result  of  (t  —  49)+  and  the  final  regression  model  is 
constructed  with  the  new  data. 

In  the  current  case,  the  filter  determines  that  the  points  at  time  49  and  55  are  not  valid  based  on 
the  change  in  slope.  These  are  not  added  to  the  expected  data  and  the  algorithm  continues  the 
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path-finding  procedure  using  the  expected  locations  in  the  x  and  y  directions.  At  this  point  of 
the  procedure,  we  know  the  target  is  traveling  from  the  departure  location  to  the  arrival  location 
on  a  straight  path  at  a  constant  speed.  Recall  that  the  optimization  model  requires  a  constant 
AOU  so  we  proceed  to  determine  the  AOU  throughout  the  target’s  path. 

The  associated  plot  of  the  AOU  and  the  resultant  MARS  model  fit  is  shown  in  Fiugre  4.19.  The 
overall  change  in  area  during  the  leg  is  less  than  the  selected  threshold  of  25%,  thus  the  AOU 
is  chosen  as  the  maximum  area  during  the  path.  In  this  case,  we  capture  90%  or  more  of  the 
probability  mass  during  the  entire  journey. 


Piecewise  Constant  AOU 


Time  (Hours) 

Figure  4.19:  Resulting  AOU  throughout  one  target’s  straight  path  and  associated  MARS  model 
fit  to  the  greedy  area  data. 


With  the  entire  track  examined  for  speed,  heading,  and  area  changes,  we  proceed  to  build  the 
final  ordinary  least  squares  model  using  the  expected  location  and  time  data.  For  this  particular 
case  with  no  waypoints  and  no  changes  in  the  AOU,  the  final  model  takes  the  following  form 


jc  =  290.5-0.205? 


(4.16) 
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and 


y  =  4.3  +  0.142?. 


(4.17) 


Based  on  the  times  in  Table  4.4  and  Equations  4.16  and  4.17,  we  determine  the  remaining 
parameters  utilizing  the  equations  described  in  Section  4.4.  Table  4.5  compares  the  results  of 
the  path-finding  model  with  the  associated  inputs  to  the  probability  model.  The  path-finding 
model  determines  the  target’s  departure  location,  speed,  and  associated  departure  and  arrival 
times  extremely  well.  When  comparing  the  destinations,  remember  that  the  path-finding  model 
uses  the  probability  data  for  the  time  periods  when  the  probabilities  sum  to  one.  Thus,  the  model 
determines  an  arrival  location  that  is  actually  on  the  water  vice  on  the  coast  as  depicted  by  the 
actual  location.  An  accurate  estimate  of  the  origin  is  required  for  the  optimization  model,  but 
an  accurate  estimate  of  the  destination  is  not  required  as  long  as  its  along  the  same  track. 


Parameter 

Actual 

Model  Result 

Speed 

15  knots 

14.981  knots 

Origin 

76.3°  W  9°  N 

76.3°  W  9°  N 

Destination 

83.5°  W  14°  N 

82.657°  W  13.414°  N 

AOU 

7200  nm2 

7582.65  nm2 

^ earliestdep 

29  hours 

29  hours 

^ min 

37  hours 

37  hours 

^ max 

64  hours 

64  hours 

f 

±4  hours 

±4  hours 

T 

33  hours 

33  hours 

Table  4.5:  Comparison  of  actual  parameters  and  path- finding  model  results  for  the  simple  target 
scenario. 


Figure  4.20  shows  the  final  model  fit  to  plots  of  the  expected  track  through  the  AOR  and  over 
time. 


4.6  Multiple  Leg  Model  and  Illustration 

Recall  the  target  with  two  legs  from  Figure  4.1.  Table  4.6  lists  the  associated  parameters  that 
we  enter  into  the  enhanced  path-finding  model. 

Again  we  show  in  Figure  4.21  several  hours  of  the  probability  mass  as  it  moves  through  the 
AOR.  The  probability  model  generates  a  heat  map  for  120  hours  of  the  target’s  journey  and  the 
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Final  Model  Fit  to  Expected  Location  for 
One  Target  on  a  Single  Leg 


Case  1 


Figure  4.20:  Output  of  the  path-finding  model,  showing  the  target’s  track  through  the  AOR  and 
the  target’s  track  in  both  the  x  and  y  direction  with  respect  to  time. 


resulting  probability  files  feed  into  the  path-finding  model.  Based  on  the  probability  data,  we 
determine  the  initial  and  final  departure  and  arrival  times  which  are  listed  in  Table  4.7. 
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Approximate  payload 

Speed 

Origin 

Departure  uncertainty 
Waypoint 

Waypoint  uncertainty 
Destination 
Arrival  uncertainty 
Departure  time 
Departure  time  uncertainty 
Intel  certainty 


2,500  kg 
15  knots 
77°  W  5°  N 
±30  nm.  Uniform 
86°  W  3°  N 
±60  nm 

92.2°  W  14.5°  N 
±50  nm.  Uniform 
10  hours 

±3  hours,  Uniform 
50% 


Table  4.6:  Parameters  for  a  multiple  leg  path  from  origin  to  destination.  Uncertainty  data  and 
associated  distributions  are  listed  for  departure  location,  arrival  location,  and  departure  time. 


Time 

Actual 

Model  Result 

^ earliestdep 

7  hours 

7  hours 

First  Tmm 

13  hours 

13  hours 

Last  xmax 

96  hours 

91  hours 

late  starr 

NA 

107  hours 

f 

6  hours 

6  hours 

T 

10  hours 

10  hours 

Table  4.7:  Case  time  information  from  path- finding  model  as  compared  to  the  actual  inputs  for 
a  target  traveling  along  a  path  with  multiple  legs. 


Next,  we  proceed  with  the  application  of  MARS  between  time  13  and  time  91.  The  resulting 
form  of  the  MARS  model  is 


xt  =  274.36  —  0.123(f  —  45)+  ±0.249(45-?)  +  (4.18) 


and 


yt  =  2.77  ±  0.052(r  -  37)+  ±  0. 170(?  -  49)  +  ±  0.058(49  - ?)+.  (4.19) 

In  this  case,  MARS  introduces  waypoints  at  hours  37,  45,  and  49.  We  examine  the  validity  of 
each  waypoint  by  looking  at  the  change  in  slope  and  the  maximum  time  for  the  entire  probability 
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-0.0 

i 


Figure  4.21:  A  one  target,  multiple  leg  scenario  produced  in  the  enhanced  probability  model. 
Snapshots  are  from  various  times  during  the  target’s  route. 


mass  to  travel  through  the  respective  waypoint.  This  time  is  then  compared  to  the  amount  of 
time  until  the  next  waypoint.  If  the  maximum  time  is  greater  than  the  time  to  the  next  waypoint, 
then  we  discard  the  next  waypoint  and  assume  it  is  a  artifact  of  the  probability  mass  moving 
through  a  turn. 

Here,  the  change  in  slope  is  significant,  thus  solidifying  that  we  have  at  least  one  valid  waypoint. 
Next  we  determine  which  of  the  three  are  valid.  We  assume  that  the  first  waypoint,  chronologi¬ 
cally,  is  valid  and  check  if  the  maximum  time  to  move  through  the  turn  is  less  than  the  time  to 
the  next  waypoint.  Figure  4.22  illustrates  the  probability  mass  at  hour  45.  The  current  waypoint 
(at  hour  45)  and  the  next  waypoint  at  hour  49  are  highlighted.  We  desire  to  determine  the  length 
of  time  for  this  probability  mass  to  move  through  the  hour  45  waypoint. 
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Probability  Mass  at  Hour  45 


Figure  4.22:  Depiction  of  the  probability  mass  in  the  one  target  multiple  leg  case.  The  prob¬ 
ability  mass  is  shown  at  hour  45,  which  is  the  time  of  the  current  waypoint.  We  rotate  this 
probability  mass  to  align  with  the  target’s  course  en  route  to  the  next  potential  waypoint. 


Specifically,  we  rotate  the  probability  mass  along  the  target’s  path  to  determine  the  maximum 
distance  between  x  coordinates.  This  rotation  aligns  the  probability  mass  with  the  target’s  track, 
allowing  us  to  use  the  target’s  overall  speed  in  the  time  calculation.  This  distance  is  then  divided 
by  target  speed  to  produce  the  time  to  complete  the  turn.  At  hour  37,  the  time  to  complete  the 
turn  is  6.8  hours,  which  is  less  than  the  time  to  the  next  waypoint  (8  hours),  thus  the  waypoint 
at  hour  45  is  kept.  The  turn  time  at  hour  45  is  7.4  hours  and  is  greater  than  the  time  to  the  next 
waypoint  (4  hours),  therefore  we  discard  hour  49  as  a  waypoint  and  conclude  the  target  path 
consists  of  three  legs  based  on  speed  and  heading  changes.  Figure  4.23  illustrates  the  rotated 
probability  mass  and  the  associated  horizontal  distance.  The  final  path  based  on  constant  speed 
and  heading  legs  contains  two  valid  waypoints  with  the  first  at  hour  37  and  the  second  at  hour 
45. 

The  AOU  analysis  follows  which  determines  if  the  path  requires  additional  legs  based  on  sig¬ 
nificant  changes  in  the  AOU.  The  associated  plot  of  the  AOU  and  the  resultant  MARS  model 
fit  is  shown  in  Figure  4.24.  In  this  case,  we  analyze  each  leg  separately  using  MARS. 

The  resultant  MARS  models  fit  three  different  linear  piecewise  segments  to  the  area  data.  The 
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Rotated  Probability  Mass  at  Hour  45 
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Figure  4.23:  Depiction  of  the  probability  mass  in  the  one  target  multiple  leg  case.  The  rotated 
probability  mass  is  shown  at  hour  45,  which  is  the  time  of  the  current  waypoint.  The  rotated 
probability  mass  aligns  with  the  target’s  course,  which  allows  us  to  calculate  the  time  for  the 
target  to  move  the  maximum  x  distance  using  overall  target  speed. 


first  segment  begins  at  hour  13  and  ends  at  hour  37.  During  this  time,  the  area  changes  by  more 
than  25%,  thus  the  model  divides  the  constant  speed  heading  leg  into  two  legs  based  on  the 
change  in  area.  The  area  on  the  other  two  legs,  from  37  to  45  and  from  45  to  91,  changes  by 
less  than  25%  and  are  not  further  divided.  Figure  4.25  depicts  the  associated  times  and  areas  of 
the  new  legs. 

With  the  addition  of  new  legs  based  on  area  changes,  the  path  now  consists  of  four  legs.  The 
new  leg  does  not  change  the  speed  or  heading  of  the  target  but  only  represents  the  location  that 
the  AOU  increases  by  more  than  25%.  Asa  result,  the  final  model  accounts  for  the  speed  and 
heading  changes  identified  by  MARS  but  it  is  unnecessary  to  insert  the  points  of  area  change. 
The  final  path-finding  model  is  of  the  form 

=  285.4  -  0.242f  -  0.024(f  -  37)+  +  0. 145  (f  -  45)+  (4.20) 


and 
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MARS  Model  Fit  to  Greedy  AOU 


Figure  4.24:  Resulting  AOU  throughout  one  target’s  multiple  leg  path  and  associated  model  fit 
to  the  greedy  area  data. 


Piecewise  Constant  AOU 


Figure  4.25:  Resulting  piecewise  constant  area  for  one  target  multiple  leg  case. 
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yt  =  5.47  -  0.050?  -  0.055(?  -  37)+  +  0.320(?  - 45)+. 


(4.21) 


The  mechanics  to  determine  the  remaining  case  parameters  are  the  same  as  those  in  Section  4.5 
and  Table  4.8  lists  the  model  output  and  the  actual  parameters.  Figure  4.26  shows  the  final 
regression  model  fit  to  the  expected  track  over  time. 


Parameters 

Leg  1 

Leg  2 

Actual 

Model 

Actual 

Model 

Speed 

15  knots 

14.843  knots 

NA 

14.843  knots 

Origin 

77°  W  5°  N 

77.015°  W  4.964°  N 

NA 

80.649°  W  4.211° 

N 

Destination 

86°  W  3°  N 

80.649°  W  4.211°  N 

NA 

83.555°  W  3.608° 

N 

AOU 

5400 

5697 

NA 

7422 

^ min 

13 

13 

NA 

25 

^ max 

46.88 

25 

NA 

37 

f 

6 

6 

NA 

6 

T 

10 

10 

NA 

25 

Leg  3 

Leg  4 

Parameters 

Actual 

Model 

Actual 

Model 

Speed 

NA 

17.175  knots 

15  knots 

14.754  knots 

Origin 

NA 

83.555°  W  3.608°  N 

86°  W  3°  N 

85.684°  W  2.763° 

N 

Destination 

NA 

85.684°  W  2.763°  N 

92.2°  W  14.5°  N 

91.264°  W  12.603c 

’N 

AOU 

NA 

9184 

10800 

10077 

^ min 

NA 

37 

46.878 

45 

^ max 

NA 

45 

96.138 

91 

f 

NA 

6 

6 

6 

T 

NA 

37 

46.878 

45 

Table  4.8:  Comparison  of  actual  parameters  and  path-finding  model  results  for  one  target  mul¬ 
tiple  leg  scenario. 


4.7  Results 

The  core  effort  of  this  thesis  concerns  the  need  to  provide  target  parameters  to  the  optimization 
model  in  Section  2.3.  Thus  far,  we  demonstrate  the  capability  to  take  intelligence  data  concern¬ 
ing  departure  and  arrival  locations,  speed,  and  departure  time  to  produce  a  probabilistic  heat 
map  over  time.  With  the  heat  map  as  input,  we  show  the  ability  to  determine  the  necessary 
target  parameters  to  build  a  search  vehicle  allocation  plan.  This  section  briefly  describes  the 
results  of  the  optimization  model  with  our  input  parameters  verse  the  results  using  the  actual 
target  parameters  for  several  target  scenarios. 
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Final  Model  Fit  to  Expected  Location  for 
One  Target  on  a  Multiple  Leg  Path 


Case  6 


Figure  4.26:  Output  of  the  path-finding  model  for  a  multiple  leg  scenario,  showing  the  target’s 
track  through  the  AOR  and  the  target’s  track  in  both  the  x  and  y  directions  with  respect  to  time. 


4.7.1  Optimization  Results 

Consider  a  scenario  with  two  targets,  where  target  1  is  an  SPSS  in  the  eastern  Pacific  and  target 
2  is  a  go-fast  in  the  Caribbean  Sea.  Target  1  follows  a  track  consisting  of  two  legs  and  target  two 
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is  on  a  track  with  a  single  leg.  Figure  4.27  depicts  the  actual  track  and  the  track  approximated 
by  the  path-finding  model. 


Figure  4.27:  Case  scenario  with  two  targets  where  target  1  is  a  SPSS  operating  in  the  eastern 
Pacific  and  target  2  is  a  go-fast  operating  in  the  Caribbean  Sea. 


Figure  4.28  depicts  the  optimal  search  plans  based  on  the  actual  track  and  the  results  of  our 
path-finding  model.  In  both  cases  the  searcher  leaves  from  home  base  and  travels  to  search  for 
target  2  and  then  to  search  for  target  1.  The  searcher’s  track  follows  the  points  in  sequential 
order,  starting  and  ending  at  home.  In  each  instance  the  search  locations  and  times  of  the 
search  are  nearly  identical,  in  addition  to  the  expected  amount  of  drugs  found.  The  fact  that 
the  two  search  plans  are  so  similar  suggests  that  our  path-finding  model  adequately  determines 
the  target’s  track  since  we  are  comparing  optimal  search  plans  with  the  parameters  input  into 
the  enhanced  probability  model  and  the  output  from  our  path-finding  model.  This  example 
illustrates  a  key  link  between  our  probability  and  path-finding  models  with  the  search  vehicle 
optimization  model. 
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Figure  4.28:  The  optimal  search  plan  based  on  the  results  of  the  path-finding  model  and  the 
actual  track.  The  two  search  plans  are  nearly  identical.  The  searcher  leaves  from  home  base 
and  travels  to  search  for  target  2  and  the  to  search  for  target  1 .  The  searcher’s  track  follows  the 
points  in  sequential  order  starting  and  ending  at  home.  The  green  track  represents  the  search 
plan  based  on  the  true  intelligence  inputs,  and  the  track  in  orange  represents  the  search  plan 
based  on  our  path-finding  model. 
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CHAPTER  5: 
Conclusion 


Illicit  drug-trafficking  is  a  major  concern  of  the  United  States  and  continued  expansion  of  drug¬ 
trafficking  significantly  affects  the  national  security  interests  of  the  United  States.  This  thesis 
addresses  a  particular  aspect  of  illegal  drug-trafficking  in  the  eastern  Pacific  and  Caribbean  Sea. 
We  focus  on  taking  intelligence  information  concerning  a  drug  trafficker’s  transit  route  to  create 
a  spatial  temporal  probability  map  and  subsequently  determine  the  trafficker’s  most  likely  route 
through  the  region.  Last,  we  address  the  accuracy  of  our  results  by  comparing  the  optimal 
search  plan  with  our  path-finding  model  as  input  and  with  the  true  case  parameters  as  input  to 
the  optimization  model. 

We  begin  with  a  fundamental  probability  model  to  account  for  uncertainty  in  the  departure  and 
arrival  locations  of  a  single  drug  trafficker  who  travels  along  a  linear  path.  Next  we  enhance 
this  model  to  account  for  additional  uncertainty  in  the  trafficker’s  speed  and  departure  time. 
Additionally,  we  build  in  functionality  to  handle  multiple  drug  traffickers,  multiple  courses  and 
speeds,  and  the  ability  to  output  results  showing  the  expected  amount  of  drugs. 

With  uncertainty  in  the  target’s  transit  route  represented  by  several  random  variables,  we  create 
a  probability  map  of  the  target’s  track  over  time.  Using  these  probabilities,  we  apply  regression 
techniques  to  determine  the  characteristics  of  the  target’s  track  needed  to  run  a  search  vehicle 
optimization  model.  The  optimization  model  requires  inputs  for  each  segment  of  the  target’s 
path  where  on  each  segment  the  target  is  traveling  along  a  constant  heading,  with  a  constant 
speed,  and  a  constant  AOU.  Departure  time  and  arrival  time  on  each  segment  identifies  where 
one  segment  ends  and  the  next  begins.  Segment  changes  occur  based  on  changes  in  speed, 
heading,  or  AOU.  Together  these  route  characteristics  represent  the  inputs  to  the  search  vehicle 
optimization  model.  Our  procedure  results  in  an  optimal  search  plan  sufficiently  close  to  one 
developed  using  the  probability  model  inputs. 

We  illustrate  the  models’  functionality  with  two  case  studies  of  varying  complexity.  First,  a 
simple  straight  path  between  the  origin  and  destination  shows  the  path-finding  model’s  ability  to 
identify  the  target’s  starting  location,  destination,  speed,  and  associated  travel  times.  Combining 
these  route  characteristics,  we  are  able  to  reproduce  the  target’s  likely  track  from  start  to  finish. 
Supplying  our  path-finding  results  into  the  search  vehicle  allocation  model  shows  similar  search 
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plan  results.  Next,  we  demonstrate  a  more  complex  scenario  with  a  target  traveling  to  the 
destination  with  a  course  change  during  the  journey.  Our  model  accurately  identifies  this  course 
change  and  reproduces  the  target’s  track.  Again,  the  resulting  search  plans  are  quite  similar. 
These  two  case  studies  are  representative  of  actual  drug  trafficker  routes  but  additional  analysis 
is  needed  to  validate  the  model’s  ability. 

Our  probability  model  and  path-finding  model  establish  additional  areas  of  research  to  aid 
JIATF  South  in  their  mission.  We  address  the  issue  surrounding  the  AOU  over  time  with  three 
different  methods  where  the  chosen  method  involves  a  greedy  algorithm.  While  this  method  re¬ 
sults  in  areas  much  smaller  than  the  other  two  approaches,  our  overall  approach  was  relatively 
informal  and  further  research  could  apply  rigorous  mathematical  methods  to  the  AOU  problem. 

We  examine  several  cases  with  various  target  characteristics  and  compare  our  path-finding 
model  results  to  the  original  probability  model  inputs.  Our  limited  comparisons  demonstrate 
promising  results  but  a  more  formal  comparison  is  necessary  to  uncover  any  modeling  shortfalls 
that  may  be  present.  Additionally,  we  compare  a  few  search  vehicle  plans  using  our  path-finding 
model  and  the  original  case  inputs.  Again,  a  more  thorough  investigation  should  be  undertaken, 
possibly  involving  design  of  experiments  or  other  factors  to  determine  where  the  model  might 
not  perform  as  well. 

We  do  not  address  cases  with  a  large  departure  uncertainty.  These  are  the  cases  referred  to  as 
"channel"  cases  in  this  thesis.  They  consist  of  scenarios  where  a  target  has  a  large  departure 
window,  on  the  order  of  ±24  hours  or  more.  The  probability  model  creates  a  long  channel 
of  probability  that  grows  from  the  origin  and  reaches  the  destination,  effectively  connecting 
the  two  locations.  In  other  words,  the  probability  mass  is  never  completely  on  the  water  and 
the  target  can  be  anywhere  along  the  channel  from  the  earliest  departure  time  to  the  latest 
arrival  time.  The  optimization  model  handles  these  cases  in  a  different  manner  and  requires 
different  inputs,  thus  follow-on  analysis  to  develop  different  algorithms  considering  these  cases 
is  necessary. 

Our  probability  model  does  not  account  for  the  impact  of  weather  on  the  drug  trafficker’s  move¬ 
ments.  A  comparison  of  the  path-finding  model  results  utilizing  input  from  our  model  and  one 
with  weather  characteristics  would  highlight  whether  or  not  those  items  are  important  factors 
for  consideration. 
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